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ABSTRACT 

This  thesis  offers  a  new  and  fundamental  approach  to  the  study  of 
ship  motion  in  a  random  seaway.   The  central  idea  is  that  the  randomness 
of  ship  motion  is  caused  by  the  forces  and  moments  of  wave  action  in  the 
six  degrees  of  freedom  of  the  ship.   These  six  forces  constitute  a  multi- 
dimensional random  process,    which  can  be  analyzed  with  recently  developed 
mathematical  theory.    The  advantages  of  this  approach  include:  (1)  The  entire 
statistics  of  the  sea  as  they  affect  ship  motion  can  be  determined  from 
measurements  made  aboard  ship,    and  (2)  Theory  is  available  for  using  the 
results  of  the  analysis  of  these  forces  to  develop  optimum  continuous  motion 
predictions  (as  would  be  useful,    for  example,    in  an  automatic  landing  system 
for  aircraft  carriers)  as  well  as  optimum  control  systems  for  the  simulta- 
neous control  of  many  stabilization  power  actuators,    such  as  fins  or  rudders. 

In  the  thesis,    mathematical  models  of  the  ship  are  presented  in  the 
Laplace  transform  and  in  the  matrix  differential  equation  forms.   The  pro- 
blems associated  with  correlation  function  and  power  spectral  density  appro- 
ximations are  discussed  and  contemporary  methods  analyzed.   An  original 
design  procedure  is  developed  which  will  allow  simultaneous  maximum -effort 
control  of  n  .stabilization  devices,   minimizing  a  generalized  quadratic  error 
criterion,    and  which  will  take  account  of  (1)  all  statistical  interrelations 
between  wave  forces,    and  (2)  the  inter-coupled  nature  of  ship  motion. 

Finally,    a  comprehensive  research  program  is  proposed  in  the  areas 
of  ship  behavior  at  sea  and  ship  motion  control  which  emphasizes  the  use 
of  multi -dimensional  random  process  theory  as  a  basic  tool. 
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CHAPTER  I. 
A  NEW  APPROACH  TO  SHIP  MOTION 

1.  1      Introduction 

This  thesis  is  concerned  with  a  new  approach  to  the  understanding 
of  the  motions  of  ships  in  a  seaway.   A  unified  theory  will  be  established 
which  is  heavily  dependent  on  the  methods  used  by  the  control  systems 
engineer,    who  has  spearheaded  much  of  the  contemporary  advances  in 
the  study  of  complex  dynamical  systems. 

Complexity  is  the  necessary  penalty  for  an  adequate  analysis  of 
the  interdependent  motions  of  a  ship  in  six  degrees  of  freedom.   With  the 
advent  of  high-speed,   large  capacity  digital  computing  equipment,    com- 
plexity and  dimensionality  alone  are  not  formidable  deterrents  to  under- 
standing,  providing  that  the  analyst  has  convenient  methods  for  expressing 
the  relations  involved  and  has  a  well-defined  objective  for  analysis. 

This  work  is  to  be  considered  as  a  practical  application  of  general 
theory  developed  in  a  companion  volume   „   the  latter  being  motivated  by 
the  dearth  of  theoretical  aids  to  understanding  of  the  random  nature  of  total 
ship  motion.   Thus  a  major  goal  of  this  thesis  is  to  interpret  and  to  illus- 
trate the  practical  application  of  the  mathematical  concepts  presented  in 
Ref.    1. 

The  hoped-for  result  of  this  work  is  to  outline  a  new  approach  to 
ship  motion  which  will  yield  dividends  in  more  complete  knowledge  of  the 
properties  of  the  ship  and  of  the  sea,    and  what  is  perhaps  just  as  important, 
will  provide  rational  design  procedures  for  control  systems  which  tend  to 
constrain  or  control  ship  motion. 

1.  2     Current  state-of-the-art 

A  very  useful  summary  of  contemporary  knowledge  on  the  subject 
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of  ship  motion  is  found  in  the  monumental  work  of  Korvin-Kroukovsky  , 

who  has  made  a  critical  evaluation  of  most  of  the  facets  of  past  research 
on  the  sea  and  its  effect  on  ships.  It  would  be  pointless  to  attempt  to  parrot 
this  work  to  any  extent,   but  there  are  many  results  of  recent  research 
which  tend  to  support  the  approach  to  be  presented  in  this  volume. 

(a)  The  hypothesis  of  linearity 

First  of  all,   consider  the  ship  at  sea,   subject  to  hydrostatic  and 
hydrodynamic  forces  which  cause  complex  motions  of  the  rigid  structure. 
Conventionally,   this  motion  is  resolved  into  translations  along  three  normal 
axes  of  a  three -coordinate  system  and  rotations  about  these  axes.   As  will 
be  illustrated  in  Section  2.  2,    these  motions  can  be  considered  as  the  solu- 
tion  of  a  six -fold  set  of  intercoupled  second-order  differential  equations, 
with  the  sea  providing  the  exciting  forces.   It  is  well-known  that  most 
systems  which  occur  in  nature  can  be  approximated  in  a  sufficiently  small 
region  of  validity  by  linear  differential  equations,    regardless  of  strong 
non-linearities  in  other  regions.   The  question  which  naturally  arises  is, 
can  the  ship  be  adequately  described  by  linear  differential  equations  in  its 
normal  operating  range  of  motion  at  sea? 

Des  pite    the  prolific  investigations  of  various  non-linearities  of 
ship  motion  which  have  appeared  in  the  literature,   it  appears  that  the 
linear  approximation  is  valid  in  many  cases  for  the  actual  ship  at  sea.   To 

support  this  hypothesis,   model  investigations  made  by  Korvin-Kroukovsky 

3  4 

and  Jacobs   ,    and  later  by  Gerritsma    in  regular  waves  showed  good  agree - 

5  6 

ment  between  experiment  and  linear  theory.    Lewis   ,    Lewis  and  Numata  , 

7 
and  Lewis  and  Dalzell     conducted  towing  tank  tests  in  irregular  seas  which 

verified  the  principle  of  superposition,   the  consequence  of  linearity,   for 
the  motion  of  the  model  with  respect  to  wave  components  at  various  fre- 
quencies.  Also,   in  Section  3.  2  of  this  work,    statistical  properties  of  mo- 
tion measured  aboard  the  USS  GYATT  (DDG  1)  in  February  1959  by  Marks 

o 

and  Durkovic     will  be  shown  to  agree  nicely  with  a  linear  hypothesis. 
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As  an  aside,   it  is  indeed  fortunate  that  the  blessings  of  linearity 
can  be  attributed  to  the  ship  at  sea,   for  a  rigorous  theory  for  non-linear 
analysis  of  complex  dynamic  devices  is  in  its  infancy,    and  when  stochastic 
processes  are  allowed  as  an  excitation  force,   the  understanding  of  such  a 
problem  is  many  years  away. 

(b)    Description  of  the  sea  and  ship  motion  as  a  random  process 

The  sea  defies  a  determinate  description.  Its  ever-changing  surface 
with  waves  rushing  here  and  there,    combining,   disappearing  --  this  would 
appear  to  be  a  phenomenon  which  is  absolutely  indescribable  in  the  large. 
However,    some  progress  has  been  made  in  delineating  the  nature  of  the 
sea  with  statistical  methods. 

The  fundamental  building  block  for  contemporary  statistical  study 
of  the  sea  is  the  unidirectional  regular  wave.   If  superposition  is  assumed 
to  apply,    which  regards  the  air -water  interface  as  a  linear  transmission 
medium,   the  sea  can  be  considered  to  be  constructed  of  a  distribution  or 
spectrum  of  waves  of  various  amplitudes,    wavelengths,   phases,    and  direc- 
tions.  Since  the  sine  wave  has  the  virtue  of  orthogonality,   has  been  conven- 
tionally used  in  the  description  of  random  processes  in  other  fields,    and 
is  an  approximation  to  the  waves  which  propagate  naturally  at  sea,    it  na- 
turally is  used  as  the  fundamental  wave  in  the  superposition  procedure. 

Some  of  the  major  contributors  to  this  ever -expanding  body  of  know- 
ledge have  been  St.   Denis  and  Pierson  ,    Neumann     ,    Darbyshire     ,    Voz- 
nessensky  and  Firsoff     ,    Barber     ,    Longuet-Higgins      ,    and  Marks 

14 
Longuet-Higgins       has  made  what  is  perhaps  the  most  appealing 

derivation  of  the  description  of  the  sea.  He  shows  that  a  generalized  auto- 
correlation function  of  the  sea  can  be  expressed  by 


<f  (x,    y,    t)    -    hm  _1 \   \       N(x  ,   y  ,   t  )    N(x-rx,   y -ry,   t  +t) 

V    y     T-*ooXYT  J  J  J  .         t        . 

Aj    7i  -X-y   -T  dx     dy     dt     (1.1) 
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Here,   N  is  the  wave  height  at  a  point,   x  and  y  are  normal  position 
coordinates,    and  t  is  a  time  variable. 

Unfortunately,   there  is  no  known  theory  which  would  permit  mani- 
pulation of  the  above  expression  or  its  possible  transform  mates  (spectra) 
—  provided  that  they  could  even  be  measured  ~~  in  order  to  reproduce  use- 
ful results  in  the  study  and  control  of  ship  motion. 

10 

Neumann       has  had  a  considerable  influence  in  the  theoretical  and 

practical  development  of  sea  spectral  theory.  By  combining  a  shrewd 

analysis  of  early  experimental  data  with  the  dynamics  of  waves  at  sea,   he 

developed  an  expression  for  the  spectrum  (or  frequency  content)  of  the 

wave -height  at  a  point 

K.  K2 


$  <i0)    =       —  e    '    w2 


-W6 

where  the  K.  terms  indicate  constants. 
1 

9 

St.   Denis  and  Pierson    hypothesized  a  directional  sea  spectra  of 

K3  -    K2 

N(a>,  9)    =     —  e        — |     cosue 

-u>6  w* 

where  ©  represents  the  direction  of  wave  travel  as  measured  from  a 
reference  direction  of  maximum  wave  travel. 

They  proposed  that  if  measurements  were  made  on  a  ship  model 
to  determine  the  response  of  the  ship,   both  amplitude  and  phase,   to  a 
si:nusoidal  wave  excitation  which  varied  over  a  range  of  frequencies  and 
a  range  of  oblique  directions  to  the  direction  of  motion,   the  response  of 
the  ship  could  be  calculated  by  weighting  each  response  function  by  the 
amount  of  excitation  at  that  frequency  and  from  that  direction,    and  summ- 
ing the  total.   This  idea  has  had  a  mjaor  effect  in  shaping  current  ideas 
about  ship  response  in  irregular  waves. 

An  alternate  approach  to  the  problem  of  defining  wave  spectra  at 

12 
a  point  was  presented  by  Voznessensky  and  Firsoff     .   They  made  ex- 
perimental measurements  of  wave  auto -correlation  functions,    which  are 
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the  inverse  Laplace  transforms  of  wave  spectra,    and  found  that  they  could 
all  be  approximated  by  an  exponentially  damped  cosine  wave,    which  will 
be  later  shown  to  agree  nicely  with  a  linear  hypothesis. 

Finally,,   it  is  necessary  to  consider  how  the  random  process  of 
wave  motion  arose  in  the  first  place.   A  number  of  investigations  have  been 

made  into  the  mechanism  of  energy  transfer  from  wind  gusts  to  wave  motion. 

2 

According  to  Korvin-Kroukovsky  ,   this  problem  is  essentially  unsolved. 

From  basic  considerations,   however,   the  wave  action  that  a  ship  sees  is 
generally  a  result  of  random  wind  action  over  great  areas  of  the  earth's 
surface  and  over  long  periods  of  time. 

(c)      Summary 

It  is  useful  to  summarize  at  this  point  the  major  facts  from  contemp- 
orary theory  which  will  be  of  use  in  the  presentation  to  follow. 

(1)  The  ship  is  in  general  suitably  described  by  the  assumption  of 
linearity. 

(2)  The  equations  most  used  for  "theoretical"  wave  spectra  are 
actually  based  on  an  analysis  of  early  experimental  data. 

(3)  The  usual  approach  in  expressing  motion  of  a  ship  in  a  random 
seaway  is  to  consider  it  as  the  effect  of  a  wave  spectra  exciting  a  linear 
system  whose  continuum  of  transfer  functions  can  be  determined  by  ex- 
periment from  ship  model  response  to  oblique  waves. 

(4)  Correlation  functions  of  wave  records  show  a  tendency  to  appro- 
ximate an  exponentially  decaying  cosine  wave. 

(5)  The  actual  cause  of  the  waves  observed  by  a  ship  are  a  very 
large  number  of  essentially  independent  wind  gusts,   the  effect  of  each 
making  a  very  small  contribution  to  the  wave  height  observed  at  a  single 
point. 
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1.  3     Outline  of  a  new  approach  to  the  ship  motion  problem 

The  heart  of  the  contemporary  approach  to  ship  motion  in  irregular 
seas  is  presented  in  Figure  1.1.  Here  the  wind  is  shown  as  a  disturbing 
variable  acting  on  the  sea  to  create  a  wave  spectrum.   The  wave  spectrum, 
in  turn,    acts  on  a  ship  to  produce  motion.   The  original  cause,   the  wind, 
is  such  a  random  phenomenon  and  is  distributed  over  such  a  large  area 
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Fig.    1.  1      Usual  approach  to  ship  motion  analysis 

that  any  attempt  to  measure  its  effect  would  be  ludicrous.   Therefore,   in 
the  large,   the  system  or  mechanism  which  converts  the  total  wind  disturb- 
ance  to  a  sea  spectrum  is  forever  unknowable.   The  spectrum  itself  is 

capable  of  measurement,   but  the  only  full  scale  attempt  to  measure  it  at 

16 
sea,   Project  SWOP     ,    was  a  very  drawn  out,    complicated,    and  expensive 

operation  to  obtain  the  statistical  properties  of  a  small  patch  of  ocean  at  a 

certain  instant  of  time.   The  continuum  of  transfer  functions  which  relate 

the  sinusoidal  response  of  the  ship  to  an  oblique  wave  can  be  measured 

in  the  towing  tank,   but  this  would  involve  a  very  elaborate  testing  program 

because  of  two  variables  involved,    direction  and  frequency.  It  appears 

that  the  study  of  ship  motion,   by  focussing  attention  on  the  causal  action 

of  a  sea  spectrum  which  is  in  practice  unmeasurabie,   has  painted  itself 

into  a  corner.   This  corner  is  the  stagnation  which  must  necessarily  result 

in  a  theory  which  depends  on  parameters  which  are  so  difficult  to  establish. 

Suppose  that  a  return  to  the  most  fundamental  principals  is  made. 
The  motion  of  a  ship  at  sea  is  mathematically  approximated  by  a  system 
of  constant-coefficient  differential  equations.  In  these  equations  appear  a 
forcing  function  of  six  dimensions  --  the  forces  along  and  moments  about 
the  three  coordinate  axes.   These  forces  represent  the  effects  of  the  random 
seaway  to  the  ship.   Fig.    1.  2  shows  the  sequence  of  actions  considered. 
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Fig.    1.2      A  new  approach  to  ship  motion  analysis 

Let  us  assume,  for  the  moment,  that  the  equation  coefficients  of 
a  particular  ship  are  known.  This  had  been  the  goal  of  a  great  deal  of 
theoretical  and  experimental  work  in  the  past  decade,   and  at  this  writing, 
it  appears  that  a  reasonable  approximation  could  be  made  on  the  basis  of 
theoretical  analysis  and  model  experiments  for  a  given  ship.  Knowing  these 
coefficients,   it  would  be  entirely  within  the  scope  of  current  technology  to 
measure  the  displacement ,   velocity,   and  acceleration  variables  of  the 
various  degrees  of  freedom  of  the  ship,    weigh  these  numerically  by  their 
respective  coefficients  in  the  differential  equations,    and  compute  the  effec- 
tive forces  on  the  ship.   That  is,    we  would  continuously  produce  on  a  ship 
at  sea  a  record  of  the  approximate  six  force  variables  which  would  be 
composed  of  suitable  weighted  values  of  measured  ship  motion.   This  six- 
fold record  would  be  a  random  process,    a  so-called  multi -dimensional 
random  process.   Providing  that  our  approximations  to  the  equations  of 
the  ship  are  valid,   we  would  have  a  record  of  the  total  random  variation 
of  the  sea  as  it  has  affected  the  ship  over  the  interval  of  time  considered  . 
It  is  important  to  note  that  the  ship  "system"  considered  here  includes 
the  hydrostatic  and  hydrodynamic  forces  which  would  act  on  the  ship  if  it 
were  suddenly  released  from  an  out-of-equilibrium  condition  in  an  other- 
wise smooth  sea.   The  random  force  inputs,   then,    are  the  integrated  super- 
imposed effects  of  the  wave  system  acting  on  the  ship. 

It  is  contended  therefore,    and  this  is  the  most  important  statement 
in  this  work,   that  the  fundamental  goal  of  the  analysis  of  ship  motion  at 
sea  should  be  to  investigate  the  statistical  nature  of  these  wave  forces  on 
the  ship. 

The  logical  question  for  the  reader  to  ask  at  this  point  is  --  How  does 
one  analyze  a  record  of  six  simultaneous  inter-related  random  processes 
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to  extract  useful  information  from  it?  A  theoretical  answer  to  this  question 
was  presented  by  the  author  in  a  companion  document     which  was  essen- 
tially motivated  by  the  particular  problem  of  analyzing  ship  motion  at  sea. 

Some  of  the  advantages  to  be  expected  from  pursuing  this  line  of 
investigation  are: 

(1)  A  catalogue  of  sea  descriptions  can  be  built  up  for  a  ship  class. 

(2)  All  statistical  data  used  can  be  collected  aboard  ship. 

(3)  A  quantitative  means  for  determining  the  effect  of  hull  coefficient 
changes  on  ship  motion  in  typical  seas  will  result. 

(4)  A  means  of  continuously  making  an  optimum  short -time  predic- 
tion  of  any  or  all  components  of  future  motion  is  provided. 

(5)  Based  on  this  approach,   a  design  procedure  is  possible  which 
yields  near-optimum  control  of  steering  and  stabilization  devices,    allowing 
(1)  for  saturation  of  the  force  transducers,   (2)  for  the  intercoupling  of  mo- 
tions in  the  ship,   (3)  for  the  statistical  relations  among  the  various  wave- 
produced  forces  and  moments,    and  (4)  for  the  possibility  of  noise -corrupted 
information  from  the  motion  sensors. 

The  following  chapters  will  develop  and  elaborate  on  the  theory  and 
practical  use  of  this  new  approach  to  the  understanding  and  control  of  ship- 
board motion. 


CHAPTER  II . 
MATHEMATICAL  MODELS  FOR  THE  SEA  AND  SHIP 

2.  1    Mathematical  tools  required 

The  study  of  such  a  complex  phenomenon  as  ship  motion  at  sea 
is  not  suitably  undertaken  with  horse -and-buggy  mathematical  methods.   In 
fact,   the  success  of  such  an  analysis  will  be  highly  dependent  on  methods 
which  can  convey  large  ideas  with  a  minimum  of  notational  difficulty. 

The  first  major  problem  which  must  be  resolved  is  the  most  efficient 
handling  of  problems  involving  the  dynamic  behavior  of  linear  systems. 
Experience  in  the  fields  of  control  and  communications  engineering  has 
shown  that  the  most  useful  expression  for  linear  systems  analysis  is  the 
transfer  function  in  terms  of  the  Laplace  transform  variable,    s.   Its  use 
avoids  the  unpleasant  task  of  convolution  when  the  response  of  a  system  is 
to  be  determined  from  knowledge  of  its  excitation.   Suppose  a  transfer  func- 
tion G(s)  of  a  system  is  given  by 

Z(s) 


G(s)    = 


P(s) 


where  Z(s)  and  P(s)  are  polynomials  in  s.  G(s)  is  the  ratio  of  the  Laplace 
transform  of  the  output  to  that  of  the  input.   The  easiest  way  of  visualizing 

the  meaning  of  this  expression  is  to  consider  the  differential  equation  which 

dn 
effectively  relates  the  input  x  and  output  y  of  the  system.   If  -r-^   is  replaced 

by  the  operator  s   ,   then  this  differential  equation  will  have  the  form 

P(s)    y    =       Z(s)    x 

Formally,   the  definition  of  the  Laplace  transform  G(s)  of  a  function 
of  time  g(t)  is  given  by 

G(s)    =  JgU)    e  "St  dt  (2.1) 

o 

Throughout  the  literature  of  naval  architecture  and,   in  particular, 
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ship  motion,   the  "jw"  notation  is  used  in  place  of  the  fundamental  transform 
variable,    s,   to  designate  the  transfer  function  of  a  physical  system.   Super- 
ficially,  this  would  seem  to  be  a  difference  of  small  importance,   but  in 
practice  it  has  demonstrated  the  extreme  preoccupation  of  contemporary 
theorists  with  the  ideas  of  simusoidal  response  of  systems,   a  point  of 
view  which  it  is  hoped  will  suffer  in  the  treatment  to  follow. 

The  second  major  mathematical  tool  which  must  be  brought  into 
play  is  matrix  theory,   often  called  the  algebra  of  higher  mathematics.   The 
use  of  matrices  permits  handling,    at  least  conceptually,  dynamic  linear 
problems  of  arbitrary  complexity.   With  six  degrees  of  freedom  of  the 
ship,   matrix  methods  are  a  necessity. 

It  is  assumed  in  the  following  that  the  reader  is  familiar  with  the 
basic  elements  of  Laplace  transform  and  matrix  theory. 

2.  2     Equations  of  motion  of  the  ship 

In  this  section  the  equations  which  govern  the  rigid-body  motion  of 
a  ship  in  its  six  degrees  of  freedom  will  be  reviewed,    and  manipulated 
into  a  form  most  suitable  for  analysis.   The  conventional  derivations  of 

the  equations  of  motion  will  not  be  presented  here,    as  they  are  available 

17 
elsewhere     ,   but  some  motivation  for  their  appearance  will  be  given. 

By  equating  the  vector  force  F  on  a  ship  to  the  time  rate  of  change 
of  momentum, 


F    =     -4r        (  m    V  )    =    m 

dt  i  w 


Idt   '  VIJ    Ivl 


and  the  vector  moment  Q  to  the  time  rate  of  change  of  angular  momentum 

a  general  set  of  equations  results  which  gives  the  three-dimensional  motion 
of  a  free  body  under  the  action  of  external  forces.  Here,   V  is  the  instantaneous 
velocity  of  the  ship,  Si  is  the  instantaneous  angular  velocity,    and  Ul_  is  the 
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angular  momentum  vector. 

A  suitable  right-handed  orthogonal  coordinate  system  is  selected: 

(1)  Under  equilibrium  conditions  at  constant  speed  U  ,  x  is  fixed 
in  the  ship  from  the  center  of  gravity  in  the  direction  of  forward  motion, 
y  is  to  starboard,    and  z  is  down. 

(2)  The  angular  displacements  -  -  <p  (roll),  ©■  (pitch),    and  ^(yaw)  -- 
are  measured  in  a  clockwise  sense  from  the  equilibrium  positions  about 
the  x,   y,    and  z  axes,    respectively. 

(3)  Linear  velocities  along  and  angular  velocities  about  the  x,   y, 
and  z  axes  are  denoted  by  u,   v,    w  and  p,    q,    r,    respectively. 

With  this  choice  of  coordinates,    six  equations  result  which  relate 
the  force  or  moment  in  a  given  direction  to  the  resulting  motion: 

»• 

F       =    m  (x  +  qw  -  vr) 

x  n 

F       =    m  (y  +  rUL   -    pw) 

...  . 

F       =    m  (z  +  pv  -  qu) 
z 

M       =    I     p  +  qr(I     -  I   ) 
x  x  z        y 

M       =1     q  +  rp(I         I   ) 
y  y  x        z 

M       =    I     r  +  pq(I     -  I   ) 
z  z  y        x 

These  equations  are  linearized  about  the  equilibrium  point,   yielding 

F       =    m  x 
x 

F       =    my  +  mUr 

y  J  o 

•  • 

F       =mz-mU     q 
z  o 

M      =    I    p 
x  x 

M       =    I     q 

y        y 
M     =   I    r 

z  z 

where  U     is  the  mean  forward  speed  of  the  ship.   The  other  terms  involving 
the  product  of  two  velocities  disappear  because  of  their  zero  equilibrium 
values. 

Next,   the  hydrodynamic  and  hydrostatic  forces  exerted  by  the  mo 
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tion  of  the  ship  are  examined.  Initially,   on  a  body  of  arbitrary  shape  it 
can  be  assumed  that  accelerations,   velocities,    and  displacements  in  or 
about  any  axis  cause  forces  in  or  about  any  other  axes.   This  would  indicate 
that  108  coefficients  in  general  would  be  needed  to  specify  the  system,    a 
staggering  thought.   A  considerable  simplification  results  in  considering 
the  real  ship: 

First  of  all,    linear  displacements  in  the  x  or  y  directions  could 
not  be  expected  to  result  in  a  force  on  the  ship.   Next,   the  symmetry  of 
the  usual  ship  form  about  its  centerline  plane  prevents  small  motion  in 
the  vertical  plane  (heave,   pitch,    and  surge)  from  inducing  forces  in  other 
degrees  of  freedom.   Also,   the  surge  or  motion  in  the  direction  of  mean 
forward  velocity  can  sensibly  be  regarded  as  being  decoupled  from  other 
motions. 

The  linear  differential  equations  which  result  from  these  considera- 
tions are  certainly  approximations  to  the  actual  case  --  the  hydrostatic 
and  hydrodynamic  forces  are  in  general  non-linear,   the  coordinate  system 
used  is  valid  only  in  a  certain  range  near  equilibrium,    and  the  assumption 
that  the  wave  excitation  forces  occur  independently  of  the  motion  of  the 
ship  is  certainly  debatable.   With  all  these  misgivings,    we  still  plunge 
ahead  and  trust  that  the  resulting  theory  will  provide  (1)  useful  results 
for  a  ship  at  sea  and  (2)  a  foundation  for  a  more  complete  theory  as  in- 
adequacies are  revealed  by  future  experimental  results. 

The  equations  of  motion  are  conveniently  expressed  in  Laplace 
transform  notation,    using  matrices, 

H(s)     v<s>]         =  Fc(s)]        +        FuMs)]  (22) 

Here  H(s)  is  a  6x6  matrix  which  expresses  the  dynamic  forces 
which  act  on  the  ship  as  a  result  of  its  motion.   V(sjis  a  six-dimensional 
vector  of  the  displacements  of  the  ship  from  equilibrium.   F  (s)Jand  FUJ(s)J 
indicate  the  forces  of  controlled  elements  and  waves,    respectively,    in  each 
of  the  degrees  of  freedom. 
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To  structure  these  equations,   the  following  components  make  up 
the  displacement  vector  v.   They  are  grouped  in  the  order  given  in  order 
to  emphasize  the  interacting  motions. 

v1  --Roll,  </>   ;  v2  --Yaw,  IfJ    ;  vg  --Sway,   y;  v^  --  Heave,    z;  v^  — 

Pitch,  G    ;  and  v.     -Surge,   x. 
b 

H(s)  is  constructed  from  three  sets  of  intercoupling  motions  -- 

H  (s),    roll-yaw-sway;  H  (s),   heave  pitch;  and  H  (s),    surge. 


H(s)    = 


t 


o    i  Hi  c* )   j    o 
6~1"~o~~,rHTi>J" 


(2.3) 


Hjts)    = 


H  (s)    - 


sV^  f   SKir 


s*(Kp  tlx)  ♦  s  KP  +  K<t>  s^r  *  s^ 

(2.4) 


f    Sl-(2^  +  a*-)    +  5  ^  + 


HJs)     =       [  S^4  V^)    +    s  Xu  J 


(2.5) 


(2.6) 


The  components  of  forces  or  moments  in  F     and  F     follow  the 

c  w 

order  of  V  above. 

This  set  of  equations  is  believed  to  be  one  possible  approximate 
representation  of  the  linear  model  for  the  ship,   but  additions  or  simpli- 
fications  will  not  affect  the  general  validity  of  the  theory  to  follow  in  future 

sections.   The  coefficients  in  the  elements  are  expressed  in  standard  nota- 

17 
tion     ,    with  the  important  exception  of  a  reversal  in  sign.   It  is  strongly 

felt  that  general  coefficients  of  motion  used  in  differential  equations 

should  be  positive,    reserving  the  negative  sign  as  an  indicator  of  such 

particular  phenomenon  as  negative  damping  or  a  negative  spring  constant, 

which  are  usual  sign  -posts  of  instability. 
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Solving  Eq.   2.  2  for  the  motion  of  the  ship, 
V(s)|   =      H^s)     /  Fc  <s)|+  F^  (s)] 


where 


H"1(s)    = 


^(s) 


0 
0 


H2"1(s) ) 


o    ;  h3  (s> 


(2,7) 


(2.8) 


Figure  2.  1  shows  the  ship  as  a  "system",   being  excited  by  wave 
and  control  forces. 

+ 


SHiP 


Fu) 


WAvJE   R)(2C£ 


COMTOOL 

force 


H  ts) 


V 


SHlP    MoTioN 


Fig.   2.  2     General  model  for  multi   dimensional  ship  motion 

Besides  the  transform  approach,   there  is  an  alternate  way  of 
characterizing  linear  systems  which  is  appealing  in  its  generality  and 

compactness  of  notation.   This  is  by  means  of  the  so-called  matrix  differ™ 

18 
ential  equation,    as  found  in,   for  example,   Bellman      .   In  this  method,   the 

primary  quantity  of  interest  is  not  the  input  or  output,   but  rather  the  state 
variables  of  the  system.   These  are  a  set  of  quantities  which,    if  specified 
as  initial  conditions,    would  be  sufficient  to  determine  tne  entire  free    re- 
sponse of  the  system  as  it  settled  to  equilibrium. 

To  motivate  this  approach,    consider  the  simple  first  order  differ- 
ential equation 

-£-     x  (t)    =    a  x  (t)    +    d    c(t) 

dt 

where  a  and  d  are  constants,    and  c(t)  is  an  independent  variable.   If  c(t) 
is  zero,    and  x(0)  is  specified,   the  solution  of  the  equation  is 


x(t)    =    x(0)    e 


a  t 


(t  >  o) 
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Now,    suppose  we  write  the  same  equation  in  matrix  notation. 


dt 


x(t) 


x(t) 


+     D     c(t) 


(2.9) 


where  A  and  D  are  matrices  with  numerical  elements,  and  x(t)  and  c(t) 
are  state  and  control  vectors,  respectively.  If  c(t)  is  zero,  and  x(0)  is 
known,   the  solution  of  this  equation  is 

(t>o)  (2.10) 


x  (t)      =    e  x(o) 


At 


from  the  theory  of  matrix  differential  equations,   e  is  an  nxn  matrix, 

known  as  the  matrix  exponential,    which  is  defined  by 

.  i 


At 


Z       A 


The  ability  to  express  the  initial  condition  response  in  such  a  simple 
form  for  linear  systems  of  arbitrary  complexity  is  a  considerable  ad- 
vantage to  the  analyst. 

The  ship  equations  of  motion  will  now  be  cast  into  the  matrix  differ- 
ential form,    and  the  numerical  elements  of  the  matrices  will  be  identified 
as  functions  of  the  ship  coefficients.   Also,   the  state  variables  of  the  ship 
will  be  defined.   In  the  following  sections  the  concept  of  state  and  initial 
condition  decay  will  play  a  prominent  part. 


From  Eq.    2.  2, 
H(s)         V(s)] 


;(si] 


+ 


(s)j 


H(s)  can  be  expressed  as  a  sum  of  three  matrices 

2 
H(s)    =      s      M    +    sB    +    R 


(2.2) 


(2.11) 


where 


M  = 


K     +  I 
P        x 

Kr 

N. 
P 

N.  +  I 

Y. 
P 

Yr 

0 

0 

0 

0 

0 

0 

K- 

v 

N. 
v 

Y.  +  m 

v 

0 

0 
0 


ZA  +  m 


M 


0) 


z. 

q 
m.  +  I 
q     y 
o     x 


0 
0 
0 
0 
0 


(2.12) 
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J 


K 


K 


K 


0 


B  = 


N 


N 


N 


v 


R    = 


Y 
P 

+m  U 
I             o 

Y 
r 

Y 

V 

0 

0 

0 

0 

0 

0 

0 

0 

0 

K4> 

0 

0 

0 

N^ 

0 

0 

yy 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

Z        Z    -m  U       0 

w        q  o 


M 


w 


M 


0 


(2.13) 


d 


(2.14) 


M         Mw  0 

z  G 


Rearranging  Eq.   2.  11  and  substituting  the  expression  for  H(s)  into 
Eq.    2.2,        s2M  V(s)]  =     ~  sB    V(sj]  -    R  V(s)]  +    F   (s)]  +    Fw  (sf\ 

Multiplying  each  side  by  M 


s2V(s)j   =     -s  M_1  B  V(sj|-  M    X  R  V(sj\+  M   1  F  (s)j+  M "1  Fw  (sj]    (2.15) 

The  following  state  variables  are  defined: 
x     =  Roll  velocity,   p 


x  =  Yaw  velocity,    r 

x  =  Sway  velocity,   v 

x  =  Heave  velocity,    w 

x  =  Pitch  velocity,    q 


x  =  Surge  velocity,   u 

x  =  Roll  displacement,  <& 

x  =  Yaw  displacement,  (f 

x  =  Heave  displacement,    z 

x  =   Pitch  displacement,  9 


The  following  identifications  can  be  made  immediately. 
2  V(sj\     =         s  X.] 
V(s)]      =       X.] 

Considering  the  product  RV(s)|   it  is  seen  that  this  is  equivalent  to 


i  =  1,    2,    ...    6 
i  =  1,    2,    ...    6 
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R  V(s)J  = 


K  (j, 

0 

0 

0 

0 

N  ^ 

0 

0 

0 

Yip 

0 

0 

0 

0 

z 

z 

Ze 

0 

0 

M 
z 

Me 

0 

0 

0 

0 

V. 


V, 


V 


=  1  ^ 


5J 


j  ■  7,    8,   9,    10 


Also,   by  definition, 


Xc 
X 


10 


X. 

x, 

X 
X, 


Assembling  these  equations  into  a  single  matrix  equation,    withj) 
indicating  a  matrix  of  all  zero  elements, 

r  s  ^ 


dt 


x 


m"1^ 


10  0  0  0  0 

0  10  0  0  0 

0  0  0  10  0 

0  0  0  0  10 


-M  =  1K 


M 


-1 


F        + 
c 


M 


U 


or,   in  general  form, 


—    x=       Ax     +     DF        +DWFW 
dt  c      c  "         w 


(2.16) 


(2.17) 


with  the  numerical  matrices  A,   D  ,   and  D^  being  defined  above  as  a 
function  of  the  coefficients  of  the  equations  of  motion. 

To  summarize,   the  ship  has  been  described  in  terms  of  a  Laplace 
transform  matrix,   H(s),    and  a  numerical  matrix  A-~  the  two  basic  des- 
criptions of  linear  systems  --in  terms  of  the  physical  parameters  and 
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hydro  dynamic  coefficients  of  ship  motion.   This  will  permit  the  efficient 
discussion  of  ship  motion  in  the  following  sections,, 

2.  3    Important  relations  in  scalar  random  process  theory 

The  purpose  of  this  section  will  be  to  present  and  analyze  some 
pertinent  relations  in  the  study  of  random  processes  and  their  effects  on 

linear  systems.   The  emphasis  here  will  be  on  viewpoint,    and  the  reader 

19 
is  referred  to  standard  texts       or  Ref.   1  for  proofs  and  extensions. 

Statistical  analyses  of  random  processes  are  made  in  order  to  de- 
tect relations  which  exist  between  signals  at  separated  instants  of  time. 
The  class  of  random  processes  considered  here  are  stationary  -     the 
statistics  do  not  vary  with  time  --  and  satisfy  the  ergodic  hypothesis, 
which  essentially  means  that  the  desired  statistics  can  be  adequately 
approximated  through  analysis  of  signals  over  a  long  interval  of  time. 
This  theory  becomes  more  applicable  as  the  process  considered  approaches 
a  Gaussian  random  process  --  which  latter  case  is  approximately  true  in 
the  ship  motion  problem. 

The  fundamental  measure  of  random  processes  is  the  correlation 
function,    ^(T),    defined  by 

y>x  <T)    =       E     /x(t)    .     y(t+r)}  (2.18) 

where  x(t)  and  y(t)  are  two  not  necessarily  distinct  random  processes, 
and  E  f .  j    indicates  the  expected  or  mean  value.   When  x  and  y  are  different, 
(-f      CO  is  known  as  a  cross -correlation  function,   and  when  y  is  identical 
to  x,     if     (T)  is  an  auto -correlation  function.   Fig.   2.  2  shows  typical  correla- 


T-* 


"H> 


Fig.  2.2  (a)  Typical  auto -correlation  function,   (b)  Typical  cross -correlation 
function. 
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Since  correlation  functions  (for  ^*>o)  have  the  appearance  of  transient 
signals  in  a  linear  system  --  that  is,    decaying  to  zero  as    ^~— =*  e=<5  it 

is  natural  to  use  their  Laplace  transform  as  an  alternate  measure  of  correla 
tion.   These  transforms  are  known  as  auto  or  cross  power  density  spectra. 

Since  correlation  functions  are  defined  for  positive  and  negative  time,   the 

20 
bilateral  or  two   sided  Laplace  transform  is  conventionally  selected 

If  jco   is  substituted  for  the  Laplace  transform  variable  s,   the  auto 
power  density  spectrum  can  have  the  connotation  of  an  amplitude  squared 
distribution  of  sinusoidal  frequency  content  in  a  signal  except  for  a  constant. 
For  a  cross  power  density  spectrum  in  jco  ,   the  information  is  conventionally 
represented  in  naval  architecture  literature  by  co  and  quadrature  spectra, 
which  are  its  real  and  imaginary  parts,    respectively. 

It  is  the  contention  of  the  present  author  that  preoccupation  with 
the  frequency  or   u>  representation  of  random  processes  tends  to  stifle 
progress  in  the  understanding  of  random  processes,    especially  when  cross 
power  density  spectra  are  considered.   For  example,   it  is  impossible  to 
acquire  much  feel  for  what  is  actually  happening  between  a  pair  of  random 
processes  from  consideration  of  their  co  and  quadrature  spectra,   yet  these 
observations  have  appeared  quite  frequently  in  recent  literature. 

As  a  more  basic  statement  of  the  worth  of  power  density  spectra 
in  the  analysis  of  random  processes,   this  thesis  will  demonstrate  that 
their  main  use  should  be  as  a  stepping-stone  in  determining  the  transfer 
functions  of  linear  systems. 

The  general  formula  which  governs  the  change  in  statistical  prop- 
erties of  signals  after  passage  through  a  linear  system  is 

<P      <s)    =         2.      5:   G.(-s)H.(s)    <£,x.y.  (s)  (2.19) 

Here,  (£  (s)  is  a  cross  power  density  spectrum  between  signals 
x  and  y  (they  may  be  identical)  or  the  Laplace  transform  of  (f  (t).  x  is 
made  up  of  n  signals,   x.,    each  passing  through  a  linear  system  with  trans - 
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fer  function  G.(s),   and  y  is  similarly  composed  of  m  y    signals,    each  of 
which  have  been  operated  on  by  H.(s).   The  set  of  input  statistics 

ffi  x.y.(s)       (  i  =  1,    2  .    .    .  n  ;     j  =  1,    2  .    .    .   m  ) 

J 

is  assumed  known  in  order  to  completely  specify  x  and  y. 

As  a  trivial  example  of  the  application  of  this  formula,    consider 
Fig.   2.  3,    where  the  input  spectra  of  x,    <f     (s),   is  specified,    and  the  out- 
put spectra  3>      (s)  and  the  cross  spectra  between  x  and  y,    ^     (s),    are 
to  be  found. 


Hi*) 


and 


Fig.   2.  3    A  simple  linear  system 
Since  x     =  x,   G     =  1  ;    y     =  x  and    H  (s)  =    H(s), 
§      (s)    =    H(-s)    H(s)     <g      (s) 

yy  xx 

1>      (s)    =    H(s)     <£     (s) 
xy  xx 


To  make  an  application  of  this  result  which  will  be  of  continuing 
significance  in  future  sections,    suppose  that  the  input  signal  x  has  an 
auto- cor relatiqn  function  which  is  a  unit  impulse  at  the  origin.   Since 
Y    (0)  =  E  I  x  (t)|,   this  impulse  would  mean  an  infinite  mean-square 
value  which  is,    of  course,   not  possible.  However  for  many  practical 
problems  the  assumption  of  an  impulse  is  valid,   just  as  in  mechanical 
problems  a  force  impulse  is  used  to  approximate  a  delivery  of  momentum 
over  a  small  finite  interval.   This  signal  is  known  as  white  noise,    and  in 
the  transfer  domain,   has  a  constant  value  at  all  values  of  oj  .   It  is  observed, 
for  example,    as  "grass"  on  a  radar  scope.   The  important  thing  about  white 
noise  is  that  is  the  ultimate  in  randomness  and  is  totally  uncor related 
with  time-shifted  versions  of  itself.   Also,   it  can  be  approximated  very 
nicely  for  use  in  analog  simulation. 
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Thus  if        $     (s)    =    1,   th 


en 


xx 


$    (s)    =    H(-s)H(s) 

yy 

It  is  known  that  power  density  spectra  of  single  random  processes 
expressed  as  rational  functions  of  s  can  always  be  factored  into  a  form 

<£(s)    =    H(-s)  H(s) 

Thus,    a  random  process  can  always  be  viewed  as  the  output  of  a 
linear  system,    denoted  by  H(s),    which  is  excited  by  white  noise.   By  postu- 
lating the  totally  random  excitation,   one  can  focus  attention  on  the  effective 
linear  system  rather  than  on  the  statistics  of  the  random  processes.   This 

viewpoint  is  not  just  an  abstraction,   for  an  hypothesis  has  been  offered 
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in  the  literature       that  most  random  phenomena  observed  can  be  traced 

back  to  an  original  cause  which  has  the  properties  of  almost-white  noise. 

An  appreciation  for  the  strength  of  this  approach  can  be  obtained 
through  considering  two  fundamental  problems  of  random  process  theory 
--the  optimum  filter  and  the  optimum  predictor.   These  new  results  were 

first  presented  in  Ref.    1,    and  illuminate  the  basic  mathematical  solutions 

22 

given  in  the  classical  work  of  Wiener 

a.      The  optimum  filter  problem 

Figure  2.4  shows  a  linear  system  W(s)  which  must  operate  on  a 
random  process  v.   The  signal  v  is  assumed  to  consist  of  signal  s  and 
noise  n,    with  the  criterion  for  optimum  performance  being  the  minimiza- 
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Fig.   2.4    The  basic  filtering  problem 
tion  of  the  mean  square  error  between  output  and  signal,    € 
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If  noise  and  signal  are  uncorrelated, 

$     (s)    =      £     (s)    +    £*     (s) 

^w  ^ss  snn 

from  Eq.   2.  19.   Let  us  assume  that  the  factors  of  the  denominators  of 
each  term  the  poles,    are  distinct.   Suppose     Q     (s)  is  factored  such  that 

<£     (s)    =    GK-s)    G(s) 

w 

with  the  poles  and  zeros  of  G(s)  all  located  in  the  left  half  of  the  complex 
plane. 


Also,   let  G(s)  be  separated  into  two  terms 
G(s)    =      S(s)    +    N(s) 


(2.20) 


where  S(s)  contains  all  the  signal  poles  and  N(s)  contains  all  the  noise 
poles,   that  is,   the  positive  poles  of  the  signal  and  noise  spectra.   The  ran- 
dom process  v  could  be  considered  to  be  constructed  from  white  noise  w 
as  shown  in  Fig.   2.  5. 
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Fig.   2.  5    Effective  generation  of  noisy  input  random  process 

Ref.    1  proves  that  the  optimum  linear  operation  on  the  random 
process  v  merely  reproduces  the  output  s  of  Fig.   2.  5  as  the  best  estimate 
of  the  actual  signal  s.   This  optimum  system  is  shown  as  a  simple  unity 
feedback  configuration  in  Fig.   2.6.   This  representation  is  also  valid  when 
signal  and  noise  are  correlated. 
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Fig.   2.6.   Simple  feedback  configuration  for  optimum  filtering 
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This  is  the  first  example  of  the  utility  of  postulating  a  model  for 
creation  of  the  random  process,   contrasting  with  the  laborious  mathe- 
matical treatment  found  in  standard  texts. 

b.      The  optimum  predictor  problem 

It  is  desired  to  produce  continuously  the  best  estimate  of  a  future 
value  of  a  random  process  in  a  mean  square  sense.   This  problem  (more 

specifically  the  prediction  of  aircraft  motion  for  fire  control  applications) 

22 

was  the  motivation  for  Wiener's    wartime  work   '    which  laid  the  foundation 

for  modern  random  process  theory. 

If  one  postulates  a  model  G(s)  for  the  process  v,    where 

G(-s)    G(s)    =      J     <s> 

then  this  transfer  function  model  is  characterized  by  a  state  just  as  any 
linear  dynamic  system  is.   These  state  variables  are  analogous  to  a  suffi- 
cient set  of  initial  conditions  needed  to  determine  the  future  output  at  v 
if  the  white  noise  excitation  were  suddenly  cut  off. 

The  optimum  predictor,    as  proved  in  Ref.    1,   is  one  which  con- 
tinuously produces  or  recovers  the  state  variables  of  the  model  G(s)  and 
weights  each  by  its  value  after  Tseconds  of  initial  condition  decay.   What 
this  means  is  that  the  best  guess  of  the  future  value  of  a  random  process 
is  made  by  neglecting  any  future  effect  of  white  noise  and  only  by  con™ 
sidering  the  future  action  of  present  signal  levels  in  the  assumed  model 
G(s),    which  summarize  the  total  effect  of  the  past  theory  of  the  random 
process  as  it  can  affect  the  future. 

The  importance  of  this  result  transcends  the  specific  application 
of  a  predictor.   Since  the  expected  value  of  future  signal  in  a  random  pro- 
cess is  found  by  allowing  the  random  process  state  variables  to  decay  as 
initial  conditions,   these  quantities  assume  prime  significance  in  the  design 
of  a  control  system  to  operate  in  the  random  environment.   These  matters 
will  be  considered  in  some  detail  in  Chapter  4. 
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2.4     Multi   dimensional  random  processes 

By  far  the  majority  of  theoretical  investigations  into  random  pro 
cesses  and  the  optimum  systems  which  operate  on  them  have  been  con- 
cerned  with  single  or  scalar  signals.   When  two  or,   in  general,   n  random 
processes  are  to  be  processed  simultaneously,   the  theoretical  aids  to 
understanding  have  been  essentially  missing.   One  of  the  major  purposes 
of  Ref.    1  was  to  meet  this  need,    especially  motivated  by  the  desire  to  examine 
random  behavior  of  ships  at  sea.   This  section  will  outline  the  problem  and 
a  general  theory  for  multi -dimensional  random  processes. 

First,   it  is  necessary  to  consider  what  a  multi-dimensional  system 
is.   We  will  concern  ourselves  with  linear  systems  which  have  inputs  and 
outputs  of  n  signals.   Fig.   2.6  shows  the  structure  to  be  considered.   Since 
each  output  q.  can  be  caused  by  any  or  all  of  the  inputs  v.,    an  entirely  gen- 
eral  representation  of  W(s)  is  the  definition 

q  .(s)    =        5     W..(s)    v.(s)        (i  =  1,   2  .   .   .  n) 

1  j*i        iJ  J 


or 


<\  (s)]  =       W(s)       V(s)] 


(2.21) 


Fig.   2.6      A  general  multi   dimensional  system. 


Each  element  of  W(s)  is  a  separate  transfer  function.   Fig.    2.  7 
demonstrates  the  configuration  of  a  2x2  matrix  system. 
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Fig.    2.  7      System  defined  by  2x2  matrix 
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Next,    we  consider  a  single  set  (or  "vector")  of  signals,   v,    and 
inquire  as  to  what  information  must  be  known  in  order  to  completely  specify 
it  as  a  random  process.  Obviously,   from  the  previous  section,   the  only 
measures  we  can  use  involve  correlation  and  thus  specification  of  the 
power  density  spectrum    _»v.v.(s)  for  i,   j  =  1,  2  .    .    .   n  is  sufficient  to 

■J 

define  the  process  statistically  for  our  purposes.   In  matrix  form,    an  nxn 
matrix  with  ijth  element  $v.v.(s),     £P     (s),    represents  v.   Similarly,   the 
cross-power  spectral  relation  ^v.  <^.(s)  is  the  ijth  element  of  a  matrix  ^      (s), 

A  general  formula  for  the  matrix  of  cross -power  density  spectra, 

(P       (s),   between  two  vector  random  processes,   x  and  y,   is 
— *Z —  £.       £=. 

35       (s)    =     __       Z     G.(-s)   <__x.x.(s)    HT(s) 


(2.22) 


where 


and 


X<s>]  ■  j5,   2ii! >   xi <s)] 

Y<sfl  =    l    HP        Y.ls)] 


The  superscript    _  j     indicates  the  matrix  transpose.   This  simple 
formula  is  analogous  to  the  scalar  case,    Eq.   2.  19. 

To  illustrate  a  simple  application  of  this  result,    consider  again 
the  system  of  Fig.   2.  6.  Here, 


and 


(B     (s)  =      W(-s)  35      <s)       WT(s)                              (2.23) 

jfjg  w 

$      (s)  =     ~§     (s)  WT(s)                                                 (2.24) 

■*-  VQ|                    vv  ______ 


Suppose  we  next  inquire  as  to  the  output  relations  of  a  white  noise 
driven  system,   since  this  concept  supplied  so  much  insight  into  the  scalar 
problem.  Here,   the  white  noise  is  an  n-dimensional  vector  of  uncorrelated 
white  noise  with  unit  power  density  spectra.    <fi    (s)  in  this  case  is  I,   the 
identity  matrix.   The  output  vector  v  of  a  matrix  system  G(s)  is,   from 
Eq.   2.23, 


5>       (s)    =    G(-s)    G    (s)  (2.25) 
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Let  us  now  ask  the  reciprocal  question  --if    $     (s)  is  found  by 
measurements  on  v,   can  a  G(s)  be  determined  which  satisfies  the  above 
relation  and  could  create  the  process  from  white  noise?  This  is  a  problem 
of  matrix  factorization  and,   once  solved,   it  serves  to  unlock  the  door 
to  the  understanding  of  multi- dimensional  random  processes. 

A  general  method  is  presented  in  Chapter  3  of  Ref.    1  which  always 
leads  to  a  solution  to  this  problem.   We  will  not  concern  ourselves  here 
with  the  details,    save  to  say  that  the  procedure  used  is  valid  without  regard 
to  order  and  consists  of  elemental  operations  such  as  polynomial  factoring 
and  solutions  of  linear  equations.   However,   the  significance  and  utility 
of  the  answer  is  of  considerable  importance. 

The  solutions  to  the  multi   dimensional  filtering  and  prediction 
problems  are  identical  to  those  cited  in  the  last  section  for  the  scalar 
case  --if  matrices  replace  the  single -dimensional  system  expressions. 
Especially  important  is  the  fact  that  the  optimum  prediction  of  future  values 
of  the  multi -dimensional  random  process  is  the  decay  of  state  variables 
determined  in  the  model  G(s). 

2.  5    Significance  of  multi -dimensional  theory  to  ship  motion  study 

A  foundation     -admittedly  sketchy  --has  now  been  laid  for  the 
presentation  of  a  new  approach  to  statistical  understanding  of  the  sea  and 
ship  motion.   The  heart  of  this  approach  is  the  visualization  of  a  random 
process,    single  or  multi -dimensional,    as  the  output  of  some  linear  system 
which  is  driven  by  white  noise.   The  randomness  of  the  random  process 
is  caused  by  the  white  noise,    and  the  memory  or  useful  part  is  caused 
by  the  linear  system. 

First,   let  us  consider  the  sea  itself.  It  is  manifestly  impractical 

1  f! 

to  monitor  a  large  area  of  ocean  for  amplitude  disturbance,   Project  SWOP 
notwithstanding.   However,    a  finite  set  of  measuring  instruments  such  as 
wave  poles  is  possible.   The  signals  from  these  devices  comprise  a  multi - 
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dimensional  random  process  and  can  be  analyzed  by  cross -power  spectra 
techniques  to  produce 
satisfies  the  equation 


techniques  to  produce  the  matrix  of  spectra,    $     (s).     The  G(s)  which 


G(-s)    GT(s)    -     3?     (s)  (2.25) 

w 

is  the  effective  system  which  carries  in  its  signal  levels  all  that  is  know- 
able  for  prediction. 

White  noise  can  be  visualized  as  a  random  sequence  of  uncorrelated 
impulses.   In  the  case  of  the  sea,   the  white  noise  is  physically  approximated 
by  the  gusts  or  random  wind  forces  which  continually  supply  energy  to  the 
sea.   The  sea  is  then  a  linear  system  which  converts  the  effect  of  a  single 
gust  into  a  small  displacement  disturbance  which  propagates  at  the  inter- 
face and  interacts  with  the  contributions  from  millions  of  other  gusts.    This 
very  gross  description  is  intended  to  make  the  point  that  white  noise  is  not 
necessarily  an  abstraction,    and  is  actually  a  major  factor  in  the  actual 
randomness  observed  at  sea. 

If  we  were  to  use  an  arbitrarily  large  number  of  wave  poles,   it  is 
possible  that  we  could  reduce  the  prediction  error  for  future  wave  motion 
at  a  given  point  to  less  than  any  desired  amount,    since  the  inclusion  of 
an  additional  signal  source  which  is  correlated  with  the  signal  in  question 
will  always  reduce  the  prediction  error. 

These  ideas  are  interesting  but  of  dubious  practical  value  since  it 
is  a  rare  physical  problem  which  is  concerned  with  instantaneous  values 
of  wave  height. 

A  second  use  of  the  multi -dimensional  theory  occurs  when  ship 
motion  is  to  be  measured.  If  the  six-dimensional  motion  of  the  ship,  q, 
is  measured,    and  the  input  forces  are  denoted  by  v,   then  from  Eq.    2.  23, 

(g      (s)    =     H   Vs)      (5       (s)     H_1'T(s) 

qq  y  w  

where  H     (s)  is  the  ship  matrix  transfer  function  as  derived  in  Section  2.2. 
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Suppose  first  that  H(s)  is  known,   from  knowledge  of  the  coefficients 

in 
given  by 


of  the  linear  differential  equations.   Then,   the  input  spectra  j)      (s)  is 


J     (s)    =    H(-s)     "§      (s)    HT(s) 

w  qq         

An  approach  which  is  less  interested  in  the  mechanism  of  ship  motion 
would  take  the  motion  statistics,    Cp      (s),    and  factor  this  such  that 

qq 

$      (s)    =    G(-s)    GT(s) 

qq  


This  system  G(s),    composed  of  a  mixture  of  force  and  motion  systems, 
has  in  its  signal  levels  or  state  variables  the  necessary  information  for 
the  optimum  prediction  of  future  motion  of  the  ship.   This  would  be  of 
considerable  practical  import  for  motion  control  devices,    and  for  perhaps 
more  exotic  uses  such  as  guided  missile  launching  installations  or  auto- 
matic landing  systems  on  aircraft  carriers. 

However  valid  this  parti  cular  approach  might  be  for  a  certain  spectrum, 
it  is  felt  that  the  relatively  constant  transfer  function  of  the  ship  should  be 
determined  and  fixed  in  any  design  procedure,    and  attention  be  paid  to  the 
real  variable  in  the  problem,   the  random  process  of  forces  on  the  ship 
due  to  wave  action. 

The  ship  coefficients  are  perhaps  best  determined  in  the  towing 
tank,   but  a  useful  check  can  be  made  for  some  of  the  major  parameters  of 
interest  --such  as  natural  frequencies  and  damping  ratios  --by  analysis 
of  the  correlation  functions  obtained  from  reduction  of  motion  data  taken 
on  the  actual  ship  at  sea.   Chapter  3  will  illustrate  an  example  of  this  par- 
ticular use  of  these  statistical  measures  in  motion  analysis  of  the  USS 
GYATT  (DDG  1). 

The  primary  emphasis,   then,    should  be  first  on  fixing  the  para- 
meters which  govern  the  differential  equations  of  motion  of  the  particular 
ship,    and  then  investigating  the  statistical  random  process  of  the  wave 
forces  on  the  ship.   The  latter  is  best  done  practically  by  simulating  aboard 


-28 


ship  the  equations  of  motion.   This  is  best  shown  by  means  of  an  example: 
Suppose  that  the  roll  moment  is  to  be  found.   From  Section  2.  2, 

MK    =    (Kp  +  Ix)  *      +    Kp  *     +  KiP  *    +  K"  V        +  Kr  *    +  Kv  y 

+  K    y    ~    M 
v  J  c 

where  <f>  ,  lb  ,   and  y  are  roll,  yaw,   and  sway  displacements,    respectively. 

These  motions  can  be  measured  quite  easily  with  accelerometers,  rate  gyros, 

etc.,    and  weighted  and  summed  as  shown.   The  control  term  M     refers  to 

the  moments  generated  about  the  roll  axis  from,   for  example,   fins  or  rudders, 

These  forces  are  probably  best  introduced  into  the  problem  from  strain 

gauges  on  the  rudder  or  fin  stock  or  other  signal  sources  which  give- an 

output  proportional  to  actual  lift  force  obtained,   not  to  the  deflection  angle 

of  the  appendage. 

The  result  of  this  measuring  and  weighting  for  all  dimensions  would 
be  a  set  of  six  random  processes  recorded  simultaneously.   The  cross - 
power  density  spectrum  3Vy(s)  would  be  calculated  between  signals,    and 

the  resulting  matrix   <2     (s)  factored  into  G(-s)  G    (s).  For  the  purposes 

°  w  

of  the  ship,   the  sea  is  mathematically  represented  over  the  interval  of 

data  taking  as  a  six-dimensional  matrix  filter  G(s)  excited  by  six  independent 

white  noise  sources. 

A  library  of  these  G(s)  representations  could  be  established.   Typical 
patterns  could  be  determined  and  useful  approximations  could  be  made  in 
their  configurations  as  various  sea  conditions  were  "tested"  by  the  ship. 
Hopefully,   the  sea  could  be  well-characterized  by  a  few  parameters  of 
the  matrix. 

A  very  important  use  of  the  G(s)  representation  is  that  it  is  in  a 
form  suitable  for  direct  instrumentation  on  an  analog  computer,    using 
commonly- available  white  noise  sources  to  excite  this  linear  system.   Thus 
the  effects  of,   for  example,    a  strong  correlation  between  pitch  and  heave 
excitation  could  be  experimentally  analyzed  for  the  purposes  of  control. 
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It  is  a  reasonable  assumption  that  G(s)  representation  of  random 
seas  would  not  differ  greatly  for  hull  forms  of  ordinary  proportions  in 
the  same  seaway.   Thus,   if  a  typical  G(s)  were  specified  after  some  ex- 
perience with  this  statistical  measure  of  the  sea,   the  effect  of  coefficient 
changes  of  a  particular  ship  could  be  evaluated  through  considering  the 
matrix 

^qq(s)    =    H_1(~s)    G(-s)    GT(s)    H~ljT(s)  (2.23) 


The  diagonal  elements  of    <j?    (s)  are  the  auto  spectra  of  roll,  pitch, 

MM 

etc.   The  mean  square  value  of  these  quantities  is  the  zero  value  of  the 
corresponding  inverse  Laplace  transform.   Thus  an  analytical  relation  is 
available  for  studying  the  magnitudes  of  motion  of  a  ship  in  a  standard 
sea  as  a  function  of  the  coefficients  in  its  differential  equations. 

With  the  increasing  use  of  digital  computers  for  complex  problems 
in  naval  architecture,   it  seems  possible  that  quite  an  ambitious  scheme 
for  hull  design  could  be  formulated.   From  a  set  of  equations  for  hull  forms 
with  several  variables,   the  coefficients  of  the  differential  equations  could 
be  calculated  by  strip  theory.   These  coefficients  could  be  used  to  evaluate 
mean-square  motion  of  the  ship  to  the  standard  seaway  G(s),    and  the  result 
used  inside  a  loop  in  the  computer  to  find  the  optimum  set  of  hull  para- 
meters. 

The  most  important  immediate  use  of  this  multi -dimensional  ran- 
dom process  appraoch  is  in  control.   This  interesting  problem  will  be  dealt 
with  in  more  detail  in  Chapter  4. 
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CHAPTER  III. 
PRACTICAL  PROCESSING  OF  RANDOM  PROCESS  STATISTICS 

3.  1    Introduction 

In  this  chapter  we  will  be  concerned  with  the  problem  of  obtaining 
useful  information  from  actual  measurements  made  on  random  processes 
in  the  form  of  correlation  functions  or  their  transform  mates,   the  power 
density  spectra.   The  theoretical  results  of  this  thesis,    and  indeed  all  ba- 
sic random  process  theory,    depend  upon  the  ability  to  approximate  corre- 
lation functions  by  exponentials  or  damped  sinusoids,   or  alternately,   to 
represent  spectra  as  meromorphic  functions  of  s,   that  is,    with  numerator 
and  denominator  which  are  polynomial  functions  of  s.   This  does  not  nec- 
essarily pose  a  significant  restriction  on  the  accuracy  of  the  theory,   how- 

19 
ever.   Laning  and  Battin     ,   for  example,   have  shown  that  correlation  func- 
tions can  be  approximated  by  the  exponential  class  of  functions  to  within 
any  desired  accuracy  level. 

In  advocating  this  approach,   the  author  hopes  to  supplant  the  well- 
known  Neumann  spectrum,   presented  in  Chapter  1,   with  a  less  rigid  model 
for  the  presentation  of  wave-height  statistics.   In  fact,   the  term  e  "**•,    which 
Neumann  used  in  approximation  of  early  experimental  spectra,    is  analyt- 
ically untractable  for  the  purposes  of  random  process  theory.   The  alter - 

12 

nate  results  of  Voznessensky  and  Firsoff  "*  however  fit  in  directly  with 

this  approach,   using  an  exponentially  damped  cosine  wave  to  approximate 
many  measured  wave -height  auto -correlation  functions. 

Also,   in  several  examples  of  correlation  functions  of  ship  motion 
to  be  presented  in  this  chapter,   the  exponentially -damped  sinusoid  will  be 
seen  as  a  fundamental  building  block.   These  particular  statistics  are  re- 
presentative of  measurements  made  aboard  the  USS  GYATT  (DDG  1)  by 

o 

Marks  and  Durcovic     in  early  1959. 
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3.  2    Correlation  function  analysis 

In  this  section,   methods  will  be  examined  for  approximating  a 
correlation  function  with  the  following  terms 


ip   <r)  = 


C.    e 

1 


-  a.f 


i  r  Ck4$  .  t+   4>.) 


allowing  the  special  case  with  (3  .  and  (p.  both  zero  to  represent  the  pure 
exponential.   Additional  but  less  frequent  allowable  terms  are  of  the  form 
'tr  e    .o     Grtsifi.T'  +   <£  .).   This  emphasis  on  decaying  sinusoidal  functions 

J  J 

is  made  because  it  corresponds  to  the  usual  observed  behavior  of  statistical 
measures  of  ship  motion  and  wave -height  records. 

The  standard  form  for  the  transfer  function  of  an  under-damped 
second  order  system  is      — - —     .  £  is  defined  as  the  damp- 


ing ratio,    and  oo 


n 


2  2 

s     +  25cj  S  +    (a) 

n  n         the  natural  or  undamped 

oct 


frequency.    A  typical  transient  of  this  system  will  contain  the  term  e        Comfit 
a   and  (3  are  related  to  the  standard  parameters,   J  and  go    ,   as  follows: 


2  2  2 

a        +    6       =     oj 
n 


and 


GO 


-  S 


n 


These  relations  are  best  visualized  with  the  ,aid  of  relations  shown 
on  a  plot  of  the  system  poles  on  the  complex  plane,    as  shown  in  Fig.    3.  1. 


Ju) 


cos-'  S 
Fig.    3.  1       s-plot    of  standard  second-order  pole 


S-PLAM& 


These  standard  parameters  have  been  adopted  because  of  their 
utility  in  plotting  frequency  response  in  control  systems  applications.   The 
damping  ratio  governs  the  rate  of  decay  of  amplitude  of  oscillation.   Fig.  3.  2 
shows  the  amplitude  ratio  of  successive  positive  and  negative  peaks  as  a 
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function  of  5  for  the  useful  ranges  of  damping  in  ship  motion. 

Another  characteristic  of  the  damped  sinusoid,   the  distortion  of 
the  normal  sinusoidal  waveform  because  of  exponential  attenuation,   is 
presented  in  Fig.   3.  3. 

When  it  is  desired  to  approximate  a  correlation  function  which  has 
been  experimentally  determined  by  one  or  more  damped  cosine  waves, 
the  procedure  becomes  quite  difficult  if  two  nearly  equal  frequencies  are 
involved,   because  of  a  "beat"  phenomenon.   In  this  case,   the  spectral  or 
transform  methods  to  be  described  in  the  next  section  become  very  attrac- 
tive.  In  the  remainder  of  this  section,    some  simple  aids  will  be  given  to 
speed  determination  of  the  characteristics  of  single  significant  sinusoidal 
components  of  correlation  functions. 

The  first  example  to  be  considered  is  an  auto -correlation  function 
tyrr  (TO  of  the  roll  angle  of  the  USS  GYATT  (DDG  1),   measured  in  state 

o 

five  and  six  seas    .   Computed  values  at  one-second  intervals  are  shown  in 
Fig.   3.4,   for  positive  values  of  7". 

First,   the  approximate  times  of  peaks  are  determined  and  averaged,   yield- 
ing 4.46  seconds  as  the  mean  time  between  adjacent  positive  and  negative 
peak  values.   Since      .  46  (3  =    n,   then  (3  =  .  705. 

The  interpolated  absolute  magnitudes  of  (fi(  T)  at  these  peak  times 
are  plotted  on  semi -log  paper.   Using  the  peak  values  eliminates  the  effect 
of  the  sinusoidal  term;  using  a  log  presentation  allows  an  exponential  to 
plot  as  a  straight  line.   Figure  3.5  shows  these  points  and  a  superimposed 
straight-line  approximation.   The  intersection  corresponding  to  the  time 

constant  of  this  line  and  the     —   ordinate  is  at  time  17.  1  seconds.   Thus 

1  e  v 

a   =    —        =  .0583.   The  standard  parameters  are  then  J  =  .082  s    and 

u     -  .707. 
n 

AT    QQ  /♦-» 

The  approximation  e  cos  .  705  7"  is  plotted  in  Fig.    3.  6,    and 

shows  quite  a  good  agreement  with  the  measured  values.  It  is  heartening 
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to  note  that  this  auto -correlation  function,    which  in  this  case  is  essentially 
a  measure  of  the  dynamic  behavior  of  the  ship  in  roll  to  a  totally  random 
excitation,    corresponds  to  a  linear  theory  even  at  large  amplitudes  of 
motion  at  sea.   The  peak  amplitude  of  roll  observed  at  the  time  these 
measurements  were  taken  was  28.5  degrees,    which  is  certainly  out  of 
the  range  of  a  small-signal  theory. 

Fig.   3.  7  presents  a  similar  approximation  to  a  pitch  auto -correla- 
tion function,    (f     (T).   The  observed  characteristics  are  5  -    .0869  and 

PP 
u>      =    .  052.   The  approximation  is  not  as  good  as  was  that  of  the  roll  func- 
tion,   which  reflects  (1)  the  non-linear  aspects  of  the  pitch  system  at  large 
amplitudes  and  (2)  an  indication  of  the  wave  excitation  statistics. 

Figures  3.  8  and  3.9  show  the  measured  cross -correlation  between 
roll  and  pitch.  Here  tf     (-T)    =  LP     if),   a  well-known  identity    .   Because 
of  the  difficulty  in  simply  expressing  these  function  in  terms  of  damped 
cosine  waves,   their  analysis  is  deferred  until  the  next  section.   However, 
it  is  important  to  note  that  the  cross -correlation  measured  between  these 
roll  and  pitch  signals  does  not  arise  principally  because  of  coupling  between 
the  roll  and  pitch  degrees  of  freedom  in  the  ship  system,   but  rather  from 
the  cross -correlation  between  the  exciting  moments  in  roll  and  pitch  from 
the  same  wave  spectrum. 

The  significance  of  correlation  function  analysis  in  the  time  domain 
is  that  dominant  frequencies  and  damping  ratios  can  be  readily  determined, 
as  well  as  an  indication  of  the  non-linearity  of  the  source  of  the  signal. 

3,  3    Analysis  of  power  density  spectra 

This  section  will  discuss  methods  of  approximating  power  density 
spectra  by  rational  transforms  in  an  efficient  fashion.   Also,    conventional 
techniques  for  digital  computation  of  spectra  from  correlation  functions 
will  be  critically  analyzed. 

The  heart  of  this  technique  is  based  on  methods  used  by  electrical 
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engineers  to  approximate  the  frequency  response  of  physical  systems  with 
transfer  functions  in  the  Laplace  variable  s.  Consider  first  the  auto  power 
density  spectrum,  &     (s),    which  can  be  represented  with  the  following 

Xa 

factors  in  numerator  and  denominator: 

K,   (+  s  +  a.)  ,  (s2  +  2  f  CO.  s  +   CO2) 

Suppose  the  logarithm  of  ^      (s)  is  taken.   In  illustration,   if 


Tjr       ,    x  K    ( -  s  +  a)     (s  +  a) 

^wlS|        •   2  2,  ,    2  2V 

(s     -2Tcj       s+oo       )(s     +  2  X  co      s+oo       ) 


XX 


n 


log  jg      (s)    =    log  K     +  log  (  -  s  +  a)  +  log  (  s  +  a  ) 


xx 

2  2  2  2 

-log  (s     -  2^  wns  +  oo      )  -  log  (s     +2JUnS+uJ 

In  general,    each  separate  pole  or  zero  has  an  additive  effect  on 
the  spectra  when  the  latter  is  expressed  in  logarithmic  form,    and  thus 
each  can  be  considered  separately. 

If  s  is  set  equal  to  j60,   the  mirror  symmetry  of  poles  and  zeros 
about  the  imaginary  axis  in  the  s=plane  for  auto  power  density  spectra 
ensures  that  each  pair  contributes  a  real  value  only.   Thus  the  magnitude 
of  each  pole  and  zero  as  a  function  of  60  is  the  item  of  interest,    and  the 
phase  angle  is  uniformly  zero  in  this  application. 

For  quick  graphical  approximation,    it  is  customary  to  plot  the 
frequency  function  on  a  logarithmic  frequency  scale.   Fig.    3.  10  shows  the 
effect  of  a  single  simple  pole, 


s  +  3    ' 


AsynPToTic 

APPROXIMATION 


Fig.    3.  10    Frequency  response  of  a  single  pole 
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Here.,   the  asymptotic  approximations  reflect  the  fact  that  for  small  CO  , 
111 


-— ^r  —7-       .   Because  of  the  logar- 

jiO  +  a/         0) 


.ZZ  —    ,    and  for  large  A)  , 
ithmic  frequency  scale,   the  high  frequency  asymptote  plots  as  a  straight 
line,    which  makes  the  entire  frequency  response  of  each  basic  building 
block  shiftable  in  frequency  without  distorting  its  shape.   As  another  example, 
Fig.   3„  11  illustrates  the  magnitude  approximation  of  a  frequency  function 


A&Hf>ToTE  S 


RfeaoEwcy 

FfWCTlOM 


(s+lr)(S+c) 


CL, 


U-  C  UJ — =? 

Fig.   3.11.   Approximation  of  frequency  function  with  poles  and  zeros 


with  one  zero  and  two  poles.  Here  the  asymptotes  are  plotted  all  at  once, 
and  the  positive  or  negative  corrections  are  added  to  approximate  the  ob- 
served frequency  function.   It  should  be  clear  that  the  use  of  what  is  eff- 
ectively a  log-log  presentation  has  resulted  in  considerable  simplification 
in  a  trial  and-error  approximation  procedure. 

The  conventional  unit  for  expressing  the  magnitude  of  the  frequency 
function,   F(j60)    ,   is  the  decibel,    defined  as  20  log10  |f  (j6))|   .   Although 
this  unit  has  a  meaning  of  power  in  audio  applications,   this  restriction  is 
not  valid  for  this  frequency  response  use  and  the  given  definition  applies 
to  any  physical  quantity  described  as  a  function  of  frequency.   For  an  asymp- 
tote of  decay  of  a  single  pole,     —    ,   the  slope  must  obviously  be  20  db  fall 
in  a  decade  of  frequency  in  00  . 

Considering  the  case  of  auto  power  density  spectra,   the  log  mag- 
nitude of  each  LHP  simple  pole  will  be  doubled  to  account  for  the  pair  of 
mirror-symmetric  poles.   A  plot  of  magnitude  of  a  simple  pair  of  poles 
of  a  spectra    is  given  in  Fig.   3.  12.   The  reminder  is  offered  that  multiplica- 
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tion  of  this  spectra  by  a  constant  K  corresponds  to  a  upward  translation 

g 
of  20  log       K  db,    and  that  replacing  s  by    —    corresponds  to  a  lateral  shift 
J.U  a 

of  the  plot  so  that  the  "breakpoint"  --or  intersection  of  the  asymptotes  -- 
occurs  at  "a"  radians/sec. 

Of  more  immediate  interest  for  the  analysis  of  wave  height  and 
ship  motion  spectra  are  the  poles  which  are  complex-conjugate  and  indicate 
sinusoidal  behavior  of  the  auto-correlation  function.   Fig.   3.  13  presents 
the  normalized  curves  for  a  spectrum  of  the  form 

§  .    I 

L-         n  n  n  n 

over  the  important  range  of  damping  ratios  for  these  applications.  Note 
that  the  asyptotic  slope  is  now  80  db/ decade,  indicating  the  presence  of 
four  poles. 

If  zeros  are  to  be  used  instead  of  poles,   the  same  curves  can  be 
used  with  a  reversal  in  sign. 

Before  using  these  tools  to  examine  several  sample  spectra,   it  is 
necessary  to  review  means  used  to  determine  these  particular  spectra 

from  their  associated  correlation  functions.   The  steps  to  be  outlined  are 

23 

given  in  the  authoritative  work  of  Marks       and  are  currently  used  at  the 

David  Taylor  Model  Basin  to  produce  spectra  from  digitally  computed 
correlation  functions. 


A  finite  duration  of  auto -correlation  function    (P     (7")  is  determined 

1  xx 

at  the  first  M  sample  points,    T  seconds  apart.   Its  Laplace  transform 

dr   (f    r  (T)    -l  e"JW'       =    coswr   -  j  sin  £01"  7  =  2   J  dT  </>     (T)  coswr 

is  approximated  by 

h-l 

2  T    Z_     tf      (kT)    cos   O.  kT       (i  =  1,    2  .    .    .  m) 

^  xx  1 

This  approximation  is  further  smoothed  by  the  operation 

xx        i  4     x  xx         i  -  1  2    ¥  xx       i  4      -^  xx     l+l 
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where  the  superscripts  s  and  u  indicate  smoothed  and  unsmoothed  respec- 
tively. 

An  auxiliary  purpose  of  the  analysis  of  these  spectra  will  then  be 
to  investigate  theeffectsrof  (1)  approximating  the  infinite  correlation  func- 
tion with  that  of  a  finite  length,   (2)  approximating  a  continuous  integral 
by  the  trapezoidal  rule,    and  (3)  smoothing  the  rough  spectra. 

The  roll  auto -correlation  function  of  the  USS  GYATT  (DDG  1),    given 
in  Figure  3.  6,    was  excellently  approximated  by  a  damped  cosine  wave, 
which  has  a  known  spectrum.   Fig.    3.  14  shows  the  rough  and  smoothed 
spectral  data  computed  with  the  above  scheme,    and  the  spectrum  of  the 
second-order  poles  with  go      and  j  computed  in  Section  3.  2.   An  excellent 
fit  is  obtained  near  the  resonant  frequency  with  the  rough  computed  points 
of  the  spectrum.   It  should  be  noted  at  this  critical  resonant  point  the 
smoothed  points  indicate  a  greater  damping  than  is  actually  the  case.   This 
is  the  real  penalty  for  such  a  crude  smoothing  operation  and  easily  ex- 
plained. 

Fig.   3.15  indicates  two  examples  of  this  form  of  data  smoothing. 
In  (a),   the  smoothed  points  tend  to  approach  the  desired  straight  line.   In 
(b),    a  peak  in  perfect  data  is  considerably  flattened  because  of  smoothing. 

R°°Gtt-  SMOOTH 


a) 


b) 

Fig.   3.15  Effects  of  data  smoothing 
With  the  resonant  peaks  expected  in  wave  height  and  ship  motion  spectra, 
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the  smoothing  operation  is  useful  in  producing  a  continuous  (if  not  true) 
curve  away  from  the  resonant  point  but  causes  a  considerable  loss  of  in- 
formation about  the  damping  of  the  system  under  investigation. 

The  points  in  Fig.    3.  14  labelled  "N"  indicate  negative  power  density 
in  the  rough  spectral  computation,    which  of  course,    could  not  occur  with- 
out imaginary  signals.   These  negative  points  at  low  frequencies  indicate 
the  effect  of  a  finite  (30  second)  length  of  auto -correlation  function  considered. 
At  high  frequencies,   the  approximation  of  the  integral  by  a  finite  summation 
causes  the  rough  spectral  points  to  be  alternately  too  positive  and  then  ne- 
gative. 

To  summarize,   the  representation  of  spectra  digitally  computed 
by  contemporary  methods  can  be  considered  essentially  accurate  only  in 
the  determination  of  resonant  frequency,   less  accurate  in  the  determination 
of  damping,    and  in  error  by  an  order  of  magnitude  at  two  octaves  away 
from  the  resonant  frequency. 

Fig.   3.  16  shows  the  rough  computed  spectra  of  the  pitch  signal  whose 
auto -correlation  function  was  given  in  Fig.    3.7.   A  good  agreement  with  the 
exact  transform  of  the  approximate  damped  cosine  wave  is  demonstrated, 
except  for  the  erratic  behavior  at  low  frequencies. 

The  above  examples  were  fairly  straight -forward  because  only  a 
magnitude  function  had  to  be  approximated.   In  the  case  of  a  cross -power 
density  spectra,   the  difficulty  is  increased  by  the  need  to  match  a  phase 
curve  as  well  as  one  of  magnitude.   This  is  identical  to  the  problem  in  con- 
trol systems  design,    when  a  transfer  function  is  to  be  found  for  a  system 
described  by  experimentally  determined  magnitude  and  phase  curves  of 
sinusoidal  response. 

This  procedure  is  strictly  one  of  intelligent  trial  and  error.   Several 

24 
good  references  are  Chestnut  and  Mayer     ,    who  treat  the  basic  problem, 

20 

and  Truxal     ,    who  discusses  more  refined  estimates  in  the  iterative  cycle. 

A  plot  of  magnitude  and  phase  for  a  single  pole  is  given  in  Fig.   3.17, 
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and  for  complex  poles  in  Figs.   3.  18  and  3.  19. 

In  the  case  of  cross  power  spectra,  one  desires  to  determine  the 
LHP  and  RHP  poles  and  zeros.  The  following  table  lists  the  correct  sign 
and    magnitude  of  the  logarithmic  magnitudes  and  phase  angles  of  Figs.  3.  17-19. 

LOG  MAGNITUDE  PHASE 

RHP  ZERO  - 

LHP  POLE  +  + 

RHP  ZERO  f 

RHP  POLE  + 

Table  3.  1.     Sign  of  magnitude  and  phase  terms 

If  the  cross  power  spectrum   $     (s)  is  to  be  found,   it  will  usually 
be  true  that  its  positive  or  LHP  poles  will  correspond  to  those  of  32    (s), 
and  its  negative  or  RHP  poles  will  be  those  of   (£      (s). 

The  roll-pitch  cross -correlation  functions  shown  in  Figs.   3.8  and 
3.9  have  a  digitally -computed  cross  spectra  given  in  Fig.    3.  20.   An  appro- 
ximation over  the  range  of  validity  appears  superimposed.   To  achieve  the 
proper  phase  response,   it  was  necessary  to  introduce  a  RHP  zero  and  LHP 
pole  pair  which  were  complex  and  mirror  symmetric.   Consideration  of 
Table  3.  1  indicates  that  the  magnitude  effect  of  these  terms  is  zero,   but 
that  the  phase  effect  is  additive,    giving  a  net  phase  shift  of  -360  degrees 
at  high  frequencies. 

The  analysis  of  cross -spectra  with  logarithmic  plots  of  frequency 
behavior  can  make  no  use  of  the  co-  and  quadrature  spectrum  records. 
Accordingly,   it  is  urged  that  magnitude  and  phase  plots  become  the  standard 
presentation. 

3.4    Conclusions 

(1)  A  review  has  been  made  of  the  normal  methods  used  for  the  approxima- 
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tion  of  correlation  functions  and  power  spectra  density  curves  with  terms 
usable  in  a  linear  theory.  No  distinct  preference  has  been  recorded  for 
either  method  --in  fact,   they  supplement  each  other  and  in  the  determina- 
tion of  the  dominant  frequency  and  damping  values  both  are  equally  efficient. 

(2)  The  method  currently  advocated  for  the  computation  of  spectra 
from  sampled  correlation  functions  causes  considerable  error  which  could 
be  lessened    with  a  more  sophisticated  computation  scheme. 

(3)  The  use  of  a  pole-zero  approximation  to  observed  wave-height  spectra 

appears  much  more  desirable,    especially  if  the  damping  ratio  $  remains 

12 
sensibly  constant  in  the  dominant  frequency  term     ,   than  is  the  preoccupa- 
tion with  the  empirically-derived  Neumann  spectrum. 

(4)  Statistical  analysis  of  ship  motion  does  not  generally  reveal  useful 
knowledge  of  the  input  wave  forces  and  moments  because  of  the  high  degree 
of  filtering  which  the  ship  system  performs  on  this  signal. 
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CHAPTER  IV. 
CONTROL  OF  SHIP  MOTION 

4.  1    Introduction  to  multi -dimensional  ship  control 

Let  us  consider  a  hypothetical  ship  which  is  equipped  with  individually 
movable  fins  at  both  the  bow  and  stern  and  a  rudder.  Motion  sensors  are 
available  to  measure  any  or  all  of  the  degrees  of  freedom  of  the  ship.   The 
design  sea  is  rough,   and  the  peak  wave  forces  exceed  the  maximum  forces 
from  the  hydrodynamic  appendages  because  of  the  latters'  finite  capacity 
and  breakdown  at  large  angles.  It  is  desired  to  minimize  the  average 
weighted  sum  of  total  squared  displacement  error  and  total  squared  velocity 
of  a  particular  point  aboard  ship  in  three  dimensional  motion,   perhaps  at 
a  critical  missile  launching  station,   by  simultaneously  controlling  each  of 
the  five  control  surfaces.   The  interactions  in  the  various  degrees  of  free- 
dom of  the  ship  are  to  be  accounted  for,    as  well  as  the  statistical  interrela- 
tionships of  the  seaway.  In  order  to  maximize  the  stabilization  per  dollar 
expended  for  the  appendage  area,   these  fins  and  rudders  are  to  be  driven 
full  throw. 

In  short,   this  problem  is  one  of  the  most  difficult  in  control  engineer- 
ing.  This  chapter  will  present  the  core  of  a  theory  which  will  allow  solution 
of  this  and  other  related  stabilization  problems  in  a  near-optimum  and 
easily-mechanized  fashion.   The  theory  will  be  applicable  to  surface  ships, 
submarines,   hydrofoil  craft,    "hydro -skimmers"  --in  short,   in  any  marine 
vehicle  problem  where  the  sea  exerts  the  disturbing  random  force  and  it 
is  desired  to  stabilize  with  several  control  variables. 

4.  2    Inadequacy  of  a  linear  theory 

Chapter  2  has  dealt  with  random  processes  and  the  optimum  systems 
to  operate  on  them.   It  is  logical  to  inquire  if  there  is  an  optimum  linear 
system  for  the  ship  control  problem.   In  a  related  example,    Fig.   4.  1  shows 
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an  input  random  process  v,    a  scalar  fixed  system  H(s),    and  a  controller 
KC(s)  which  acts  to  reduce  the  output  error  e.  In  this  simple  case  it  is 
obvious  that  including  an  infinite  gain  in  the  controller,   (ie:  K  — >   <:>0     ) 

5/STEM 


IT 
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Hfc) 
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><< 
—  / 
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l<  CCsj 

^ 

Fig.   4.  1    System  subject  to  disturbance 

will  cause  the  output  or  error  to  be  zero  (if  the  resulting  system  is  stable) 
since 

E(s)  H(s)  fc  1  asK-^-o 


V(s) 


1  +  K  C(s)  H(s) 


K  C  (s) 


Unfortunately,   in  the  ship  motion  problem  infinite  gain  is  impossible 

in  the  "controller",    whose  output  is  viewed  as  the  force  of  the  fins  or  rudder. 

This  infinite  gain  would  imply  an  infinite  force  output  for  a  finite  error,   which 

is  patently  not  the  case.   The  real  limitation  for  ship  motion  control  is  the 

peak  force  which  can  be  exerted  by  the  hydro -dynamic  control  surfaces. 

25 
From  Newton     ,    a  theory  exists  which  gives  the  optimum  linear 

controller  when  its  output  is  constrained  to  have  a  certain  average  power 

output.   This  appears  superficially  to  be  a  valid  approach,   but  actually  it 

is  forcing  the  physical  problem  to  fit  the  theory,    rather  than  vice-versa. 

It  further  limits  the  corrective  force  available  instead  of  making  intelligent 

use  of  the  known  limitation. 

The  designer  of  a  ship  pays  a  considerable  penalty  in  space  and 
weight  for  large  fins  and  rudders  and  their  associated  high-power  actuating 
devices.   Because  of  this  high  price,   he  naturally  desires  to  maximize  the 
stabilization  obtained  from  his  installation,    even  at  the  expense  of  using 
quite  sophisticated  signal  processing  equipment,   including  perhaps  many 
motion  sensors  and  a  digital  computer.  With  this  assumption,   we  must 


55- 


discard  our  hopes  for  a  completely  linear  theory  since  experience  has 

20 
shown       that  non-linear  control  systems  in  general  have  much  better 

performance  in  practice  than  linear  systems,   if  the  former  are  adequate- 
ly designed. 

If  an  error  or  the  difference  between  the  actual  state  of  the  ship 
and  the  desired  equilibrium  position  exists,   the  control  system  will  attempt 
to  reduce  this  error  by  exerting  corrective  effort.   It  is  intuitively  appeal- 
ing that  this  corrective  action  would  be  "better"  at  every  instant  if  its 
magnitude  were  greater.   That  is,    a  maximum -effect,   full-throw  or  "bang- 
bang"  control  system,    which  at  every  instant  uses  the  maximum  lift  force 
possible  at  every  hydrodynamic  surface,    will  be  postulated  in  the  develop- 
ment to  follow.   A  linear  mode  of  operation  in  an  arbitrarily  small  region 
of  system  error  will  be  permissible,   however.   This  full-throw  type  of 
control  has  been  shown  to  be  optimum  in  a  class  of  problem  which  is  con- 
cerned with  the  settling  of  a  free  system  to  equilibrium  in  minimum  time. 

4.  3    A  multi- dimensional  approach  for  saturating  controls 

The  design  philosophy  to  be  presented  in  this  section  was  first 
given  in  Ref.  1,  and  is  basically  an  extension  of  the  so-called  "Second 
Method  of  Lyapunov",    a  very  general  theory  of  dynamic  behavior  of  linear 

and  non-linear  systems  which  is  currently  receiving  considerable  atten- 

26 
tion  in  the  literature      .   A  summary  of  this  approach  is  presented  first, 

followed  by  a  more  complete  development  of  the  theory. 

In  general,    we  first  attack  this  problem  by  specifying  a  quadratic 
criterion  of  instantaneous  error  for  the  ship,    which  could  be,   for  example, 
the  weighted  sum  of  all  the  various  squared  displacements  and  velocities 

(ie:  state  variables)  existing  at  a  certain  time.    More  formally  we  postu- 

18     T 

late  an  instantaneous  measure  of  error  with  the  quadratic  form       x     Q  x, 

where  x  is  the  state  vector  of  the  ship  from  Section  2.  2,    and  Q  is  a  non- 
negative  definite  symmetric  matrix. 
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In  control,    we  are  interested  in  the  future  as  well  as  the  present, 
and  a  measure  of  the  future  error  of  our  instantaneous  state  is  quite 
naturally  given  by  the  integral  of  the  instantaneous  error  criterion  over 
all  future  time.   Since  a  random  process  is  exciting  the  system,   the  ex- 
pected value  of  future  components  of  x  is  used  in  place  of  a  determinate 
x„   Also,   future  values  of  x  are  a  function  of  the  random  process  state 
variables,    as  well  as  those  of  the  ship.   By  considering  y  to  be  the  state 
vector  of  a  linear  system  which  includes  both  the  ship  and  the  random 
process  generating  model,   the  Second  Method  of  Lyapunov  is  used  to  de- 
termine a  symmetric  matrix  P,    such  that 

T  |  T 

y       Py=        J       dtxQx 

T 

Thus,    at  every  instant  of  time,    a  single  number  y    Py  is  a  measure 

of  all  knowable  future  error,   neglecting  future  values  of  control  inputs. 
The  various  control  signals  are  based  on  the  requirement  that  the  negative 
rate  of  change  of  this  measure  of  future  error  be  maximized  continuously. 
The  advantages  of  this  approach  are: 

(1)  The  theory  becomes  more  valid  as  the  sea  becomes  rougher 
and  the  saturated  outputs  of  the  force  transducers  become  relatively  feeble 
with  respect  to  the  wave  excitation. 

(  2)  The  design  procedure  is  simple  to  execute  and  results  in  a 
completely  linear  system  except  for  the  saturation  constraint. 

(3)  The  resulting  system  is  guaranteed  to  be  stable. 

(4)  The  multi -dimensional  and  inter-coupled  aspects  of  the  problem 
do  not  appreciably  complicate  the  solution. 

Now,   to  consider  in  detail  the  aspects  of  this  procedure. 

First,    it  is  necessary  to  select  the  instantaneous  criterion  of  error, 

T 

x    Qx.   This  is  the  area  in  which  the  skill  and  knowledge  of  the  designer 

is  applied  with  best  effect.   As  in  all  optimization  procedures,   the  critical 
and  often  least  understood  part  of  the  design  is  deciding  what  is  really 
wanted. 
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One  rational  choice  is  for  Q  to  be  a  diagonal  matrix  with  positive 

T  <c  2 

elements,  X.  .     x    Qx  would  then  be  ^  A.  x.    .   The  X.  could  be  selected 

T     i  X      1     1  1 

so  that  x    Qx  is  approximately  proportional  to  the  total  energy  of  the  ship. 

1  2 

Since  kinetic  energy  is  represented  by    —    M  V     and  potential  energy  by 

12 

—  K  x  ,    with  K  an  effective  spring  constant,   the  X weighting  terms  could 

be  selected  proportional  to  the  acceleration  and  displacement  coefficients 

in  the  differential  equations  relating  the  velocity  and  displacement  state 

variables,    respectively,    neglecting  the  cross -coupling. 

Another  potentially  useful  selection  of  the  matrix  Q  would  represent 

error  by    X   .  X    +  u   V   .  V     ,   the  squared  disturbance  displacement  plus 
a      a     *™»     a      a 

the  weighted  value  of  squared  velocity  in  three  dimensions,   at  a  single 
point  "a"  on  the  ship.   To  obtain  this  measure  easily,    suppose  that  the 
vector  t  locates  the  point  a  from  the  center  of  gravity  of  the  ship.   If  the 
vector  vrepresents  the  instantaneous  disturbance  velocity  of  the  center 
of  gravity  -     that  is,   the  actual  velocity  less  the  mean  forward  velocity, 
the  vector  O  represents  the  angular  displacement  (for  example,   the  roll 
displacement  §  would  be  a  vector  along  x  axis  with  a  magnitude  equal  to 
(j)  in  radians),   the  vector  x  is  the  disturbance  displacement  of  the  center 
of  gravity,    and  the  vector  &J  is  the  total  angular  velocity  of  the  ship,   then 
the  motion  of  the  point  a  is  given  by 

x      =    x    +    0  x  r 
a 

\f      =    v    +     4)  x   r 
a 

for  small  angles  of  displacement. 

Let  the  state  of  the  ship  be  defined  as  in  Section  2.  2, 

x  =  Roll  velocity,   p  x      =  Surge  velocity,    u 

x  =  Yaw  velocity,    r  x_    =  Roll  displacement,    <f> 

x  =  Sway  velocity,   v  x      =  Yaw  displacement,    ^ 

x  =  Heave  velocity,    w  x      =  Heave  displacement,   ? 

x  =  Pitch  velocity,    q  xm~  P^ch  displacement,    0 

Then,   from  evaluation  of  the  vector  equations, 
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where  r  ,    r  ,    r     locate  the  point  of  the  ship  to  be  stabilized  from  the  ship 
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Having  determined,   by  any  means,    a  measure  of  instantaneous 

T 
error,   x    Qx,    it  is  now  necessary  to  evaluate  the  integral  of  this  quantity 

over  all  future  time.   To  do  this,    we  use  particular  results  from  the  Second 
Method  of  Lyapunov      .    First  of  all,    let  us  define  the  system  under  con- 
sideration.  Fig.   4.  2  illustrates  the  creation  of  the  random  forces  and 
moments  of  the  sea  from  a  white  noise  driven  model  G(s).   These  forces 
excite  the  ship,   H     (s),    and  a  controller,    sensing  the  motion,    adjusts  its 
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various  control  variables  with  full  throw. 


Section  2.  2  showed  how  the  matrix  differential  equation 


gj-       X     =     A    X 


+    D  c    F   c    +     Dw    FW 


was  formed  for  the  ship  alone.  If  the  concept  of  system  is  enlarged  to 
include  G(s),    which  of  course  has  its  own  state  variables,   then  the  differ- 
ential equation  of  this  new  system  can  be  written  as 


dt         7    =    B    y    +    D    F 


(4.1) 


where  y  now  contains,   besides  the  ship  state  variables,   the  state  variables 
of  G(s).   F  is  the  excitation  vector,   having  controlled  force  inputs  and 


white  noise  as  its  components. 
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Fig.   4.  2    The  multi -dimensional  control  problem 

If  y  has  for  its  first  ten  variables  the  corresponding  state  variables 

T 

of  x,    which  is  certainly  allowable,   the  instantaneous  error  criterion  y    Qy 

can  be  formed  by  using  the  Q  matrix  discussed  earlier  in  the  first  ten  rows 
and  columns,    and  filling  in  the  remainder  of  the  matrix  with  zeros.   This 
is  equivalent  to  saying  that  we  assign  no  instantaneous  error  to  a  state 
variable  of  the  sea,    since  we  cannot  control  it. 

If  the  relationship 


r     t  t 

J    y      Qydt    =     y      Py 
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is  to  be  satisfied,   then  by  differentiating  the  above  equation  we  obtain 


.00 


Now,   in  the  computation  of   J    y       Q  y  dt  we  are  concerned  with 
the  expected  value  of  future  y,    disregarding  future  control  action.   Hence 

d 


dt    y  -   B  y 


Substituting, 


"yQy    =    yB1Py    +    y1PBy    =    y       (B     P    +    P  B)    y 

Therefore, 

BTP    +    PB    =    -  Q 

The  elements  of  the  symmetrical  matrix  P  can  be  solved  for 

numerically  by  the     -z independent  equatibns  in  the  above  matrix 

equations,    where  n  is  the  dimension  of  y. 

The  criterion  to  be  used  to  fix  the  sign  of  the  control  variable  is 

T 

that  the  measure  of  future  error,   y       P  y,    is  to  have  the  maximum  possible 

rate  of  decrease  at  every  instant.   In  other  words, 


d 
dt 


{'***)*{*    ,T}    pm/p|=^}\ 


T 


is  to  have  a  maximum  negative  rate  of  change.   But 

4?        =         B  y    +    D  F 

dt  J 

Substituting, 
i    {yTPy}=    2      {yTBT    +    FTDT}py 

The  effect  of  the  control  variables  should  be  such  that  F 
is  as  negative  as  possible.   This  will  obviously  occur  when  F.,   the  i      con- 

trol  variable,   is  at  its  maximum,    and  has  the  opposite  sign  of  the  i      com- 

T 

ponent  of  D     P  y.   That  is, 

F      =     -      J  F    max  j      sgn    I  DT  P  y   J    .  (4.2) 
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{dT  p  y] 


where  sgn    x    =    +  1,     x   >  O 

=     -  1,    x    <  O 

Eq.   4.  2  is  the  desired  control  law. 

T 

The  numerical  coefficients  of  -  D    P  provides  a  weighting  of  the 

measured  values  of  displacement  and  velocity  in  the  various  degrees  of 
freedom  of  the  ship,    and  also  of  the  state  variables  of  the  random  process 
generating  model,   G(s). 

The  measurement  of  these  latter  state  variables  poses  an  interesting 
problem.   First  of  all,    with  the  use  of  accelerometers,    rate  gyros,    etc., 
it  is  feasible  to  measure  the  motion  of  the  ship.   Thus,   the  equations  of 
motion  of  the  ship,    summarized  by  His),    allow  direct  reproduction  of 
the  wave  forces,    F    ,   by  suitable  weighting  of  the  outputs  of  each  of  the 
motion  sensors. 

This  six-dimensional  random  process  will  have  been  previously 
analyzed  by  statistical  methods  to  obtain  a  matrix  of  power  density  spectra. 
Suppose  that  there  is  a  small  but  finite  noise  in  each  of  the  signals  re- 
presenting the  wave  forces,    which  would  certainly  be  the  case  in  practice, 
and  an  approximation  to  the  statistical  nature  of  this  noise  is  available 
(probably  a  white  noise  assumption)  . 

Even  if  there  are  steady  state  sinusoidal  components  of  noise, 

which  might  occur  with  sustained  vibration  of  the  ship's  structure,   for 

example,   this  approach  is  still  applicable  and  results  in  filters  which 

tune  out  the  offending  sinusoidal  components.   By  representing  the  input 

spectra  of  the  wave  forces   as  measured  by   <P     (s)    =    V     (s)    +    <fc?     (s), 

w  ss  nn     " 

where    $     (s)  represents  the  signal  or  force  spectra,    and    <§    (s)  the  un- 
wanted noise  spectra,   then    Q     (s)  can  be  factored  into  G(-s)  G    (s)  by 
the  methods  of  Ref.    1. 

G(s)  is  separated,   by  partial  fraction  expansion,   into  S(s)  and  N(s) 
terms  respectively,    where  for  the  assumption  of  uncorrelated  white  noise, 
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N(s)  will  be  a  numerical  matrix.   Then,   the  multi -dimensional  feedback 
system  shown  in  Figure  4.  3  is  constructed. 
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Fig.   4.3.   Optimum  multi   dimensional  filter  used  to  produce  state 
variables. 

This  shows  the  configuration  of  an  optimum  multi -dimensional 
filter,    as  was  discussed  in  Sections  2.  3  and  2.  4.   In  this  application,   the 
output  or  filtered  wave  force  signals  are  not  used,   but  the  inner  signals 
of  S(s)  are  the  desired  state  variables.   Fig.   4.4  outlines  the  proposed 


scheme  of  control. 
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Fig.   4.4.   A  full-throw  control  system  for  multi -dimensional  ship 
stabilization 

Here,  K     weights  the  outputs  of  the  motion  sensors  according  to 

— -  -1 

the  equations  of  motion  of  the  ship  to  calculate  the  external  forces.  _N_ 

and_S_act  as  an  optimum  filter  to  reproduce  the  state  variables  of  the  wave 
force  random  process.  K     weights  the  state  variables  according  to  Eq.  4.  2, 
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T 
with  the  matrix,    -  D    P. 

30 

Broniwitz       has  conducted  an  analog  computer  investigation  of  a 

particular  scalar  problem,    stabilization  of  a  second  order  system,    using 
this  theory.  His  results  show  a  superiority  of  this  approach  as  the  system 
damping  ratio  becomes  low  and  the  ratio  of  maximum  control  force  to 
r.m.s.   disturbance  force  becomes  small,    compared  with  other  "optimum" 
design  procedures.   As  this  domain  of  superiority  is  especially  applicable 
to  ship  motion,   considerable  promise  is  offered  for  the  multi-dimensional 
case. 

Thus,   to  summarize,   the  multi   dimensional  full-throw  controller 
is  made  up  of  signals  which  it  is  practical  to  measure  and  performs 
essentially  a  linear  computation,    except  for  the  binary  nature  of  the  output 
control  commands.   A  fair  amount  of  analog  measurement  and  calculation 
is  called  for,   but  the  basic  problem  itself  is  a  very  complex  one,   perhaps 
among  the  most  difficult  encountered  in  control  systems  engineering. 

4.  4    Some  problems  of  practical  multi -dimensional  ship  motion  control 

The  design  presented  in  the  previous  section  is  in  a  sense  a  central- 
ized one  --it  considers  the  ship  as  a  whole  and  makes  continuous  decisions 
as  to  how  each  control  variable  should  best  be  switched.   In  this  section 
we  will  consider  several  problems  which  arise  in  the  practical  mechaniza- 
tion of  this  approach: 

(1)  Practical  force  transducers,    such  as  hydraulically  driven  fins 
and  rudders,    require  a  finite  time  to  reverse  the  direction  of  maximum 
force  and  cannot  "switch"  instantaneously. 

(2)  The  actual  hydrodynamic  force  exerted  on  a  fin  depends  con- 
siderably on  the  relative  direction  of  motion  of  the  local  water  flow  and 
not  entirely  on  its  mechanical  displacement  with  respect  to  the  ship. 

(3)  Practical  motion  sensors  can  be  expected  to  deliver  a  signal 
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which  is  only  an  approximation  to  the  true  ship  motion  because  of  internal 
noise,    cross -coupling  effects,    and  local  vibration  levels. 

In  considering  the  first  objection,    we  assume  that  the  actuating 
equipment  is  sufficiently  powered  to  allow  full  travel  in  a  time  equal  to 
or  less  than  the  mean  switching  interval.   If  this  assumption  is  not  satis- 
fied,  then  the  system  is  rate  limited  rather  than  amplitude  limited  and 
the  design  premises  are  shaky  indeed.   In  this  case,    certainly  the  capacity 
(ie:  area)  of  the  fins  would  not  be  efficiently  used  and  either  a  reduction 
in  their  size  or  an  increase  in  the  power  rating  of  the  fin  mover  would 
appear  to  offer  a  more  economically  optimum  design  choice. 

The  output  of  the  weighting  operation <D     P  yj.,   it  will  be  recalled 
from  Section  4.3,    is  the  multiplier  of  the  i      control  force  which  determines 
the  rate  of  change  of  the  error  criterion.   When  it  has  the  value  zero,   ob- 
viously the  force  exerted  at  that  instant  contributes  nothing  to  the  reduc- 
tion of  this  criterion  and  might  as  well  be  zero.   That  is,    since  a  finite 
time  is  required  to  go  from  full  up  to  full  down,    a  near  optimum  choice 
of  when  to  commence  the  downward  travel  would  be  so  as  to  make  the 
force  pass  through  zero  at  the  same  time  as  <=  D     P  y  >■ .  does.   This  would 
minimize  the  effect  of  the  necessary  small  values  of  force  required  in 
a  finite-time  passage  through  zero  by  weighting  them  with  small  values 
of  multiplier. 

XL 

Thus,   the  signal  for  the  i      force  variable  to  commence  changing 
sign  should  be  the  prediction  of  sgn  I-  D     P  y|  .  at  T.  seconds  in  the  future, 
where  2  T.  is  the  time  required  for  full  travel  of  the  i      fin.   This  will  be 


given 


f       T  B  T      "1  BT 

by  sgn-(-  D     P  e        iy|,    where   e        i    is  the  matrix  exponential 


defined  in  Section  2.  2,    and  B  is  given  in  Eq.   4.  1.   The  switching  computer 

th  T  BT 

K_  shown  in  Fig.   4.  4  would  then  have  an  i      row  of  -  D    P  e        i  . 

The  second  problem,   that  of  relative  flow  effects,    is  not  so  easily 
handled.   The  trouble  is  pictured  in  Fig.   4.  5,    which  shows  the  lift  break- 
down of  a  typical  hydrodynamic  surface.   Suppose  that  a  feedback  signal 
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is  available  of  the  actual  lift  on  the  fin,    which  is  customary  in  contemporary 
fin  installations.   If  the  maximum  force  available  and  commanded  is  B,    and 


PIS  PLACE  MENT      RELATIVE       TO 
/NSTANTANEOUS       FLOW 


Fig.   4.5.   Ambiguous  force -displacement  relation  in  fin  breakdown 

the  actual  force  measured  is  given  by  the  dashed  line  in  Fig.   4„  5,   then 
it  is  not  clear  whether  the  relative  displacement  is  at  A  or  C,   and  the 
direction  for  corrective  motion  is  unknown  without  (for  example)  additional 
pressure  distribution  measurements. 

This  problem  is  quite  difficult,    and  typical  of  a  general  "peak -hold- 
ing" class  of  control  problem.   While  the  mechanization  of  the  solution  to 
this  interesting  dilemma  will  not  be  elaborated  upon  further,   one  possibility 
would  employ  a  small  high-frequency  sinusoidal  "dither"  on  the  fin  --or 
on  a  lift  flap  --  and  the  phase  of  the  resulting  sinusoidal  measured  force 
would  distinguish  between  the  case  of  point  A  and  C.   If  the  measured  al- 
ternating force  is  in  phase  with  the  dither,   then  fin  breakdown  has  not 
occurred,    and,    in  fact,   this  sinusoidal  force  signal  represents  a  valid 
error  signal  for  a  positioning  control  system  for  the  fin  since  it  goes  to 
zero  at  the  desired  peak  output. 

The  compensation  of  the  noisy  outputs  of  motion  sensors  is  a 
standard    problem  of  instrumentation  engineering.   If  the  noise  can  be 
statistically  approximated  by  a  power  spectra,   the  theory  in  the  previous 
section  for  optimum  recovery  of  signal  information  is  immediately  appli- 
cable.  Also,    a  convenient  help  in  this  respect  would  be  the  installation 
of  sharply  tuned  filters  at  the  instrument  outputs  which  would  reject  sinus- 
oidal components  at  dominant  hull  vibration  frequencies  excited  by  slamm- 
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ing,   for  example. 

4.  5    Adaptive  measurement  of  wave  force  spectra 

The  key  assumption  in  the  previous  theory  is  that  the  random  process 
which  represents  the  six-dimensional  wave  forces  on  the  ship  is  stationary. 
This  is,   of  course,   not  true,    since  ship  course  and  speed  changes  vary 
the  spectra  as  felt  by  the  ship  and  since  the  properties  of  the  sea  change 
over  a  period  of  hours.   However,   the  time  variation  of  the  parameters 
of  the  sea  is  very  slow,    so  that  the  process  can  be  termed  a  quasi -stationary 
random  process  and  the  previous  theory  is  valid,   providing  that  the  current 
statistics  of  the  sea  are  known. 

There  are  three  possible  approaches  to  this  problem: 

(1)  Since  the  statistics  of  the  sea  are  not  known,   neglect  them  and 
their  associated  control  theory,    and  carry  out  designs  by  more  conventional 
(if  less  effective)  means. 

(2)  Consider  the  worst  case,    statistically,    which  the  ship  could 
encounter  (that  is,    excitation  at  the  resonances  of  the  ship)  and  design 

for  this  condition,    accepting  the  degradation  of  performance  at  less  severe 
sea  states. 

(3)  Determine  the  dominant  characteristics  of  the  sea  spectra, 
and  attempt  to  measure  them  continuously  and  alter  the  control  system 
to  reflect  these  changes. 

Some  possibilities  of  the  third  approach  are  to  be  discussed  in 
this  section. 

It  is  difficult  to  discuss  intelligently  the  properties  of  the  wave- 
force  random  process  without  experimental  verification.  However,   there 
have  been  some  statistical  analyses  of  the  sea  which  indicate  considerable 
promise  in  adaptive  measurement  of  wave  forces. 

12 

Voznessensky  and  Firsoff    "  made  an  extensive  series  of  measure- 
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merits  of  wave  heights  at  sea  using  a  wave  pole,    and  performed  auto -correla- 
tion analyses.   They  found  that  the  correlation  function  was  strikingly  ap- 

-ap- 
proximated in  all  cases  by  a  damped  cosine  wave,   K     e         cos(3T  .   From 

the  discussions  of    Chapter  II,   it  is  clear  that  a  generating  model  with 

the  transfer  function 

K2 


2  2 

s     +  2  Cw'      s+w 
J      n  n 


excited  by  white  noise  --  a  standard  second  order  system  --  could  be  used 
to  reproduce  the  process.   The  important  result  of  these  measurements 
was  that  the  ratio  of  —    remained  essentially  constant  at  .  21,   despite  the 
wave -heights  and  directionality  of  the  particular  wave  system  measured. 
The  constant  ratio   /1        fi2~    t   is  the  non-dimensional  damping  ratio    5 
in  the  standard  form.   This  constancy  appears  to  be  a  natural  property  of 


waves  at  sea,    and  tremendously  simplifies  the  measurement  of  the  spectra, 
as  will  be  discussed. 

The  wave  spectrum,   if  5"  is  constant,   is  totally  specified  if  the  mean- 
square  value  of  the  signal  and  w     is  known.   The  mean-square  value  of 
the  signal  is  easily  determined  by  a  squaring  and  filtering  operation  of 
an  electrical  signal  representing  the  process. 

The   co     value  can  be  found  through  a  process  of  zero -counting,    or 

determining  the  average  number  of  zero -crossings  of  the  signal  per  unit 

27  28  2Q 

time.   From  the  work  of  Rice     ,   Ehrenfeld,    et  al     ,    and  Steinberg  '    it 

can  be  shown  that  for  a  standard  second-order  system  excited  by  white 

noise,   the  expected  number  of  zero-crossings  per  second  is     un 

The  possibilities  of  evaluating  the  dominant  characteristics  of  a 
random  process  by  squaring  and  averaging,    and  zero-counting,   is  a  power- 
ful tool  for  continuously  evaluating  the  parameters  of  the  sea.   It  is  hypoth- 
esized that  the  dominant  quantities  of  the  multi -dimensional  wave  force 
spectra,    since  they  arise  from  wave  motion,    can  be  measured  by  these 
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simple  operations.   As  a  simplified  example,   the  effective  frequency  of 
excitation  of  the  pitch -heave  system,    and  that  of  the  roll  excitation,    could 
be  measured  continuously  by  zero-counting  the  calculated  forces  as  re- 
produced in  Fig.   4.4. 

The  ability  to  make  good  statistical  estimates  of  the  parameters 
of  the  input  random  process  without  the  need  for  a  full-scale  correlation 
or  spectral  analysis  is  a  considerable  boon  to  the  control  systems  designer. 
It  is  felt  that  a  valid  approximation  to  the  third  alternative  at  the  first  of 
this  section  would  be,   on  the  basis  of  future  experimental  measurements 
of  spectra  at  sea,   to  fix  many  of  the  parameters  as  sensibly  constant  and 
to  measure  some  of  the  most  important  characteristics  --  in  particular, 
"frequencies"  of  excitation,    w      --  and  design  a  control  system  which  would 
be  nearly  optimum  in  the  sense  of  this  chapter  and  which  would  be  contin- 
uously "adapted"  to  changes  in  ship's  course  and  speed  and  sea  spectra. 

4.  6    Summary 

The  ship  motion  control  problem  is  one  of  great  complexity  because 
of  the  inter-coupling  of  ship  motion,   the  correlation  between  the  various 
wave  forces,    and  the  saturation  of  control  transducers.   This  chapter  has 
presented  one  solution  to  this  problem  which  has  the  virtues  of  near  optimal 
performance,    and  the  possibility  of  mechanization  with  a  reasonable  amount 
of  computing  equipment.  Needless  to  say,   only  the  germ  of  the  idea  appears 
here,    and  the  real  test  of  its  validity  must  await  accumulated  measurements 
of  ship  force  spectra  and  subsequent  analog  computer  simulation  of  the  con- 
trol system. 
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CHAPTER  V. 
CONCLUSIONS  AND  RECOMMENDATIONS 

5.  1    Introduction 

The  penalty  for  making  a  frontal  attack  on  such  a  complex  problem 
as  ship  motion  in  a  random  seaway  is  necessarily  a  compromise  between 
a  desire  to  examine  all  facets  in  detail  and  a  need  to  communicate  the  es- 
sential ideas  clearly.   This  chapter  will  first  attempt  to  outline  the  major 
viewpoints  which  have  comprised  the  backbone  of  this  thesis.   Then,    a 
specific  program  of  proposed  research  to  implement  these  ideas  will  be 
presented. 

5.  2    Major  ideas 

The  basic  idea  which  threads  through  all  past  chapters  is  that  the 
random  wave  forces  and  moments  on  a  ship  at  sea  are  the  dominant  quan- 
tities of  interest  in  any  study  of  the  complicated  behavior  of  multi-dimens- 
ional  ship  motion.   With  the  development  of  suitable  theory  in  Ref.    1  to 
analyze  the  statistical  interrelations  among  these  forces  and  moments, 
all  useful  information  regarding  (1)  average  motion  levels  ,   (2)  short  term 
motion  prediction,    and  (3)  optimum  control  decisions  can  be  found.   In  short, 
a  path  is  now  opened  which  will  explain  all  that  can  be  explained  about  ship 
motion  within  the  limits  of  current  mathematical  knowledge  and,    what  is 
perhaps  more  significant,   the  needed  measurements  can  be  easily  made 
with  standard  instrumentation  and  data  reduction  methods. 

Chapter  II  has  briefly  shown  that  the  ship  can  be  described  in  the 
two  basic  models  for  linear  systems  --  the  matrix  transfer  function  in 
Laplace  transform  notation  and  the  matrix  differential  equation,   the  latter 
using  the  significant  concept  of  "state".   The  idea  of  a  white  noise  driven 
model  has  been  introduced,    and  the  statistical  measures  of  power  density 
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spectra  have  been  interpreted  as  means  for  determining  this  linear  model. 

The  practical  problem  of  approximating  experimental  spectra  with 
rational  Laplace  transforms  has  been  examined  in  Chapter  III.   Besides 
presenting  some  graphical  aids  in  this  task,   it  has  been  demonstrated 
that  methods  now  in  use  for  statistical  data  reduction  using  a  digital  com- 
puter introduce  unnecessary  error  in  the  final  spectra. 

Also,   it  was  shown  that  a  considerable  degree  of  statistical  correla- 
tion exists  between  the  roll  and  pitch  motion  of  a  ship  at  sea,   the  major 
portion  of  which  is  hypothesized  to  arise  from  the  correlation  between  the 
respective  wave  excitation  forces. 

Chapter  IV  has  developed  in  some  detail  a  proposed  scheme  of 
simultaneous  control  of  many  force  transducers  to  stabilize  a  ship  at  sea. 
With  the  virtues  of  near-optimality,   full=throw  control,    case  of  instrumen- 
tation,   stability,    and  consideration  of  all  interactions  in  the  ship  motion 
and  the  input  forces,   this  scheme  offers  considerable  promise  of  extract- 
ing near  the  maximum  in  stabilization  from  fins  and  rudders. 

30 

Preliminary  analog  computer  investigation  reported  elsewhere 

gives  substantial  support  to  this  design  concept  for  ship  control.    The 
peculiar  problems  of  relative  water  velocity  effects,   hydraulic  rate  satura- 
tion,   and  adaptive  spectral  measurement  have  been  discussed. 

5.3    A  proposed  research  program 

The  following  is  a  proposed  program  of  needed  and  unified  basic 
research  in  the  areas  of  ship  behavior  at  sea  and  ship  motion  control. 

A.     Ship  parameter  measurement 

Determine,    experimentally  and  analytically,   the  significant  co- 
efficients of  the  linear  differential  equations  (as  a  function  of  forward 
velocity)  which  govern  the  total  motion  of  the  particular  ship  type  under 
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investigation.    Use  free  or  constrained  models  with  a  suggested  greater 
use  of  initial  condition  response.   Also,    statistical  analysis  of  measured 
ship  motion  at  sea  will  yield  useful  verification  of  the  small-scale  results. 

B.  Wave  force  measurement 

On  a  ship  at  sea,    or  on  a  free  model  in  a  random  wave  generating 
facility,    make  sufficient  weighted  motion  measurements  to  infer  --  through 
the  assumed  linear  differential  equations  of  motion  --  the  wave  force  or 
moment  excitation  in  each  of  the  six  degrees  of  freedom.   An  alternate 
experiment  could  make  use  of  a  completely  rigid  model  in  a  random  wave 
tank,    with  the  desired  forces  measured  directly. 

C.  Immediate  wave  force  analysis 

From  the  six-fold  record  of  wave  forces>  conduct  a  complete  cross  - 
spectral  analysis,    either  by  digital  computation  of  correlation  functions 
and  subsequent  transformation,    or  by  direct  analog  measurement.   Approxi  = 
mate  these  spectra  by  rational  transforms.   Factor,   preferably  with  a 
digital  computer,   the  matrix  of  power  density  spectra  with  the  method 
introduced  in  Ref.    1.   The  result,    a  6x6  matrix  of  transfer  functions, 
completely  characterizes  the  random  nature  of  the  sea  and  can  be  imme- 
diately used  in  analog  computer  studies. 

D.  Extended  wave  force  analysis 

Collect  a  number  of  matrix  descriptions  of  the  sea.   Attempt  to 
determine  relatively  invariant  parameters,    as  well  as  the  most  significant 
varying  ones.   Define  a  typical  or  design  sea,   possibly  on  a  worst-case 
basis  for  the  ship  under  consideration. 

E.  Relation  to  wave -height  directional  spectra 

Conduct  theoretical  and  experimental  analyses  of  the  relation 
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between  the  observed  statistical  properties  of  wave  forces  on  a  ship  or 
model  and  the  directional  spectra  which  characterizes  the  sea. 

F.     Multi   dimensional  control 


Using  the  matrix  model  developed  under  steps  C  and  D  above,    con- 
duct analog  computer  investigations  with  the  following  objectives: 

(1)  To  evaluate  the  control  theory  outlined  in  this  thesis,    as  well 
as  to  stimulate  other  approaches  to  the  multi - dimensional  stabilization 
problem. 

(2)  To  determine  the  importance  of  using  knowledge  of  the  input 
random  process  as  an  aid  to  control,    as  well  as  the  penalties  associated 
with  imperfect  or  no  knowledge  of  it. 

(3)  To  determine  the  optimum  balance  between  large  fast  power 
actuators  of  rudders  and  fins  and  large  surface  areas  of  these   appendages, 
and  also  the  reduction  in  stabilization  because  of  rate  saturation  and  lift 
breakdown. 

G.     Ship  motion  analysis 

From  accumulated  knowledge  of  wave-force  spectra,    determine 
the  significant  parameters  of  the  ship  which  affect  seaworthiness  and 
emphasize  undesirable  motion. 

H.     Adaptive  measurement  of  seas 

Conduct  analytical  and  experimental  investigations  into  adaptive 
measurements  or  "tracking"  of  significant  parameters  of  random  processes, 
with  the  goal  of  providing  a  shipboard  control  system  with  a  self-adjusting 
capability  as  sea  conditions  change. 

5.4    Conclusions 

The  above  program,    which  unifies  the  central  themes  of  this  thesis, 
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is  believed  to  offer  a  promising  and  practical  roadmap  into  the  vital  areas 
of  ship  motion  and  control.   Many  problems  and  an  unbelievable  amount 
of  data  remain  to  be  analyzed  and  digested,   but  no  valid  approach  to  such 
a  complex  problem  as  this  can  avoid  such  a  large-scope  investigation. 
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